arXiv:1502.06296vl [astro-ph.GA] 23 Feb 2015 


A radio-polarisation and rotation measure study of the Gum Nebula 

and its environment 


C. R. Purcell^^, B. M. Gaensler^, X. H. Sun^, E. Carretti^, G. BernardR’^, M. Haverkorn® 
M. J. Kesteven^, S. Poppi®, D. H. F. M. Schnitzeler^ 

and 

L. Staveley-Smith^’® 


ABSTRACT 

The Gum Nebula is 36°-wide shell-like emission nebula at a distance of only ~ 450 pc. It has been 
hypothesised to be an old supernova remnant, fossil Hii region, wind-blown bubble, or combination of 
multiple objects. Here we investigate the magneto-ionic properties of the nebula using data from recent 
surveys: radio-continuum data from the NRAO VLA and S-band Parkes All Sky Surveys, and Ha data 
from the Southern H-Alpha Sky Survey Atlas. We model the upper part of the nebula as a spherical shell of 
ionised gas expanding into the ambient medium. We perform a maximum-likelihood Markov chain Monte- 
Carlo fit to the NVSS rotation measure data, using the Ha data to constrain average electron density in 
the shell ne- Assuming a latitudinal background gradient in RM we find rig = 1.31 q '4 cm“^, angular radius 
douter = degrees, shell thickness dr = 18.5l4;4pc, ambient magnetic field strength Bq = 3.912^2 

and warm gas filling factor / = 0.3dQ4. We constrain the local, small-scale 260 pc) pitch-angle of the 
ordered Galactic magnetic field to -1-7° ^ p ^ -1-44°, which represents a significant deviation from the 
median field orientation on kiloparsec scales (~ —7.2°). The moderate compression factor X = at 

the edge of the Ha shell implies that the ‘old supernova remnant’ origin is unlikely. Our results support 
a model of the nebula as a Hii region around a wind-blown bubble. Analysis of depolarisation in 2.3 GHz 
S-PASS data is consistent with this hypothesis and our best-fitting values agree well with previous studies 
of interstellar bubbles. 

Subject headings: radio continuum: general - radio continuum: ISM - surveys - magnetic fields - techniques: 
polarimetric - ISM: individual objects (Gum Nebula) 


1. Introduction 

Observations of atomic Hi, molecular clouds and 
photo-dissociation regions in the Galaxy have shown 
that gas in a wide range of environments is gath¬ 
ered into spheres, bubbles or shell-like structures (e.g.. 
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Jackson et al.ll2005 Ghurchwell et al.l2006l:lMcGlure-Griffiths et al 


2009lj . Most of these objects are formed by physical pro¬ 
cesses associated with the evolution of high-mass stars 
(> 8M0). During their time on the main-sequence, 
such stars emit high fluxes of ultra-violet photons and 
fast winds of particles that ionise expanding H ii regions, 
evacuate low-density cavities and sweep gas into shells 
in a ‘snow-plough’ effect. At the end of their lives the 
stars eject their outer layers, before exploding as super¬ 
novae, driving strong shocks into the interstellar medium 
(ISM). OB-type stars generally form in clusters, so the 
combined action of stellar winds and coeval supernova 
explosions can gi ve rise to ‘supersh ells’, hundreds of par¬ 
secs in size (e.g.. IMoss et al. 2012ll . Supernovae and su¬ 
pershells are thought t o power the circulation of material 
into the Galactic halo (Dove et al. 200(11 : Reynolds et al 


2 OO 1 L IPidoDrvhora et alT 2007h and play a leading role 
in sculpting the fractal structure of gas in the Galactic 
disk. With energies greater than ~ 10^^ ergs, supernovae 
are also believed to be the main driver of turbulence in 
the disk ( McGrav fc Sno\^ 1979^ and have been shown 
to trigger new episodes of star-formation when shocks 
over run and compress pre-existing c lumps of molecular 
gas ( Reipurth 19831 : Oev et al. 2005 ). 




































Spheres or bubbles of plasma also present an excellent 
opportunity to probe conditions in the ISM. Studies of in¬ 
dividual objects can yield information on the conditions 
in the medium into which they are expanding and on 
their progenitors. Of particular interest are recent works 
that have used observations of Faraday rotation to derive 
the magne to-ionic properties of bubbles and their local 
ISM (e.g. Kothes fc Brown 2009t Whiting et al. 20091 : 


Harvev-Smith et al.l l201ll : ISavage et al.l 1201,11 1 . As the 
bubbles expand they interact with the ordered magnetic 
field of the Galaxy, compressing the ambient medium 
and the field parallel to the shock front. The resulting 
field geometry is a function of the pre-existing field con¬ 
figuration and the rate of expansion, leading to a unique 
rotation measure (RM) signature on the sky. While ob¬ 
servations of RMs from extra-galactic point sources yield 
only the average line-of-sight field strength, modelling 
the RM-signature of a supernova remnant or H ii region 
is one of the few ways to measure the local magnetic field 
on scales of a few hundred parsecs. Such measurements 
are essential ancho r-points for studies of t he large-scale 
Galactic field (e.g.. iKothes fc Brownll2009ll . 

If the line-of-sight magnetic field strength is known, 
Faraday rotation is a good tool to measure density jumps 
in the ionised ISM as oc RM. The density jump 
present at the shell boundary is a key indicator of the 
type of object powering the expansion. For example, 
mass and momentum conservation in the radiatively 
cooling shock-fronts of old supernova remnants (SNR) 
are expecte d to lead to very high density jumps at their 
boundaries ( Shull fc Draine 1987t l. In contrast, the ion¬ 
isation front of an evolved H ii region created by a clus¬ 
ter of B-type stars would expand at the local sound 
speed (typically ^ 10 kms“^), creating only a slight den¬ 
sity jump. The level of compression and magnetic field 
strength in the ionised gas have profound implications 
for whether star formation is triggered or suppressed by 
the passing shock. 

In this paper we present a study of one of the most 
prominent bubbles in the southern sky: the Gum Neb¬ 
ula. We start in (fl.ll bv reviewing the literature on the 
nebula, summarising its properties and theories of ori¬ 
gin. In 521 we introduce the datasets and images used in 
this work. We go on to describe our ionised shell model 
and analysis techniques in 21 We present the results 
of fitting the model to the rotation measure data in ijH 
where we derive strength and direction of the magnetic 
field, and the density jump across the edge of the nebula. 
Discussion and further analysis of depolarisation at radio 
wavelengths are presented in 21 Finally, we present our 
conclusions in 21 and suggest future avenues of investi¬ 
gation. 

1.1. The Gum Nebula 

The Gum Nebula is one of the largest o ptical emission 
nebulae in the southern sky (Gum 1953) ■ R has an ap¬ 
proximately circular morphology w ith an angular diame¬ 
ter of 36° ( Ghanot fc Sivanl[l983 1. its centre is thought 
to lie at a distance of ^500 pc from the Sun and its radius 


is ^ 130 pc ( Woermann et all 12001 ). Orig inally discov¬ 
ered in large-area photographic plates by IGu 3 (ll9,52lf . 
it dominates modern H a maps of the southern Galac¬ 
tic plane (e.g. . Dennison et al.lllO^ Gaustad et al.llioOl 
Haffner et al. 2n0,3h . 


1.1.1. The environment of the Gum Nebula 

Gonsiderable controversy exists in the literature on 
the origin and evolution of the Gum Nebula. In part, 
this is because the nebula straddles the mid-plane of the 
Galaxy and its footprint encompasses a large number 
of overlapping objects: Hu regions, supernova remnants 


vestigations of the nebula fe.g.. Gum 19561: 

Brandt et al. 

197ll 1 Alexander et al. 

1971: BeuermannI 1973: iRevnolds 

1976bl: Weaver et al.l 19771: Vallee & Bignelll 1983^ were 


limited by the paucity of observations covering the entire 


a clear pictnre ( 

3 .g.. Sahu & Sahu 1992: Sahu & Sahu 

199,31: 

Duncan et al. 19961: Revnoso & Dubned 1997: 

Woermann et al. 

2001: Stil & Tavlod 20071). Fig. fTl 


presents an ann otated H a image of the nebula in Galac¬ 
tic coordinates ( Finkbeiner 20031) illustrating the prin¬ 


cipal structures identified to-date. The upper third of 
the nebula is relatively free of confusing sources except 
at (/, 6 ) ss (268°,-1-13°) where the Ha shell overlaps 


the Antlia supernova rem nant (SNR, iMcGulloueh et al 


2002tll^obelli et al.ll2014[) . The lower two-thirds contain 
the majority of confusing objects, only some of which are 
directly associated with the Gum Nebula. 

The energy budget of the nebula is dominated by 
the output of early-type stars: f Puppis; an 04f star, 
and 7 ^ Velorum; a Wolf -Rayet star of type WG 8 
with an 07.51 companion ( De Marco fc Schmut j 1993) . 
7 ^ Velorum is embedded in the Vela OB2 association, 
which contains a furt her 81 B-type stars at a mean dis¬ 
tance of 415 ± 10pc ( de Zeeuw et al. 1999ll . The com¬ 
bined flux from C. Pup and 7 ^ Vel is capabl e of maintain¬ 


ing th e ionisation state of the Gum Nebula (|Weayer et al 


19771 1. howeyer. Vela OB2 also appears to be creating a 
smaller bubble within the Gum Nebula - the IRAS Vela 


Shell (IVS). The IVS was identified by ISahu fc Sahu 
1 1993h as a radius = 7.5° ring-like structure in the 100 /rm 
IRAS Sky Suryey Atlas centred on Vela OB2 It is asso- 
ciated with a thick shell o f Hi ( Dubner_e£^ 19921) and 
swept-up molecular gas ( Ghurchwell et al r il996ll , and 
has been i nterpreted as a wind -blown bubble driven by 
VelaOB2 ( Sahu &: Sahu 1993 ). Two further OB-stars 
are in the field: the 06 star GD-47 4551 lies well beyond 
the Gum Nebula at a distance of ^ 1300 pc, too far to 
be significantly interacting with the nebula. The star 
HD 4 9798 is a sub-dwarf 06 binary ( Bisscheroux et akl 
19971) located just outside the nebula (d = 600 ± 100 pc) 
and is observed to be emitting a wi nd that is distortin; 
the lo wer shell of the Gum Nebula ( Revnoso fc Dubne: 

1993 ). 

A strin g of Hii regions (e .g., RGW19, RGW27 & 
ROW 33; Rodgers et a.l.lll96B) are visible in Fig. [I] bi¬ 
secting the nebula along the Galactic plane. Most 



























































































































































+ 20 . 0 ' 


+ 10 . 0 °- 


0) 

^ +00.0°^ ' 

2 

u 


u 

JD 

(D 

(J3 


- 10 . 0 °- 


- 20 . 0 °- 


Antlia SNR 


280 . 0 ° 


Gum Nebula 


420 
360 

300 

240 



2 ^ 0 ° 


260 . 0 ° 

Galactic Longitude 


250 . 0 ° 


240 . 0 ° 


Fig. 1.— An Ha image of the Gum Nebula (|Finkbeineijl2003l l annotated with significant objects identified in the literature. 
OB-type stars ar e mark ed with '+’ symbols, the Vela pulsar with a square and the kinematic centre of the nebula derived by 
IWoermann et all 1 I 2 OOIII with a triangle. Boundaries of gas-disk s or shells are marked us ing coloured circles. See tll.llfor details. 
Refer ences for the an notated objects are a s follows: Antlia SN RriMc Cullough et al.| [2002l : Giant HI Ga s Disk: iReyno so fc D ubned 
I 1997 I: H i Gas Shell: I.Dubner etZlll992l: IR AS Vela She ll: ISahu fc Sahulll99l '^RGW-Hii re gions: iRodeers et al.lll960l : Vela 
Pulsar: iRadhakrishnan fc Manchester! 19691 : GD-474551: lReedll2003l : HD 49798 and C Puppis: Ivan LeeuwerJl2007l : -7^ Velorum: 
iDe Marco &: Schmut jfl9^ . 

















































of these are associated with the Vela molecular ridge 
(VMR), a concentration of molecular c louds beyond the 


Gum Nebula at a dis ta nce of 1 - 2 koc (iMav et al 


iMurphv fc Mavl ri99lh . Revnoso fc Dnbneil (jlQQ? ) dis- 


19881: 


covered a massive (1.4 x 10® M 0 ) H i gas disk correspond¬ 
ing to the optical outline of the Gum Nebula and specu¬ 
late that this may be the signature of the expanding rear 
wall of the nebula on the VMR. 

The bulk gas motions and excitation conditions of the 
Gum Nebula shell have been measured via spectroscopy 
of optical emission lines. Spectra of H a, [N i i] A6584 , 
[On] A5007 and [Her] A5876 taken bv iRevnoldsI ( 1976a ) 
suggested that much of the emitting gas is confined to 
a shell of radius = 125 pc with an expansion velocity ^ 
20kms“^, a thickness L « 15 —30pc and a temperature 
of 11300K. The expansion velocity was later updated to 
a value ~ 10 kms“^ and the excitation conditions in the 
Gum shell measured to be consistent with an H ii region 
( Wallerstein et al. 1980l : Srinivasan et al. 1987 1. 

The kinematics of the Gum Nebula have also been 
studied via observations of cometar y globules: dense ac¬ 
cretions of molecular g a s and d u st (iHawa rden fc Brand 
1976t Sandavist 1976 : Zealey 1979 : ReipurthI 1983 


Sahu et alT 1988 ; Sahu fc Sahu 19921 1993h The most 


comprehensive analysis of the nebula kinematics was 
carried out by IWoermann et al.l (| 200 l[l who found that 
the best-fitting model of the neutral gas (including OH 
masers and molecular clouds) was an asymmetric ex¬ 
panding shell whose front face is expanding faster than 
the rear (14kms“^ versus 8.5kms“^). The runaway 
0-star ( Puppis was within < 0.5° of the expansion cen¬ 
tre (Z = 261°, b = —2.5°) approximately ^ 1.5Myr 
ago, leading to speculation that its companion star ex- 
ploded, ejected C Pupp is and created the Gum Nebula. 
Woermann et al.l ( 200 ill question if the arc of H a emis¬ 


sion at 6 > 10 ° is part of the nebula, as it lies offset 
in Galactic latitude from the best-fitting neutral shell. 
However, we note that the upper part of the nebula is 
not well sampled by any of the datasets used. Only one 
data-point from that study (from a diffuse molecular 
cloud) lies at 6 > 10 °, so fits to the upper nebula are 
p oorly constrained. 

iDuncan et al.l ( 1996h estimated that synchrotron emis¬ 
sion is responsible for only 10 to 20 percent of the total- 
power from the nebula in their 2.4 GHz single-dish map, 
which covered the interior region (| 6 | < 5°). The hydro¬ 
gen radio re combination lines H156a and H139a were 
detected by IWoermann et al.l ( 2000 ll at four positions 
confirming that bremsstrahlung is the dominant radio 
emission mechanism in the upper shell. 

1.1.2. Origin of the Gum Nebula 

Four different models have been proposed in the lit¬ 
erature to explain the origin and evolution of the Gum 
Nebula: 

1. A large and moderately evolved ('^ 10® yr) Hii re¬ 
gion, i.e., a Stro mgren sphere excited by C Pu ppis 
and 7 ^Velorum ( Gumlll95^ Beuermannlll973 1. 


2. An old (> 1 Myr) supernova remnant that has 
now cooled and whose shell is subsequently be- 
i ng ionised by the e arly type stars in the interior 
( Brandt et al. 1971 : Alexander et al ][i97l|). 


3. A stellar wind bubble blown by C Puppis with 
help from 7 ^ Velorum and the Vela OB 2-association 
( Reynolds! 1976bl : IWeaver et al. 1977ll . 


4. A supershell resulting from the combination of 
multiple supernova explosions and photoionising 
effects powered by a si ngle stellar association 
( Revnoso fc PnbnCT 19^. 


Any successful model must explain the thin ionised shell 
{R/dr ^ 15), low expansion velocity (~ 10kms“^) and 
optical spectra consiste n t with low excitatio n conditions 
( Srinivasan et al.l 1987 : Sahu fc Sahul Il99^ . Glassical 
Stromgren sphere H ii regions exp and at approx imately 


the observed velocity ('^ 4kms ^ 
not produce a shell structure . A scaled version of the 


Lasker 1966h but do 


supernova model of iGhevalieil (jl974[l can produce a bub¬ 
ble of the correct size, but we would then expect to 
see significant radio synchrotron emission from the edge 
of the n ebula and this is no t detected in observations 
to date ( Haslam et al. 1982ll . The old supernova rem¬ 
nant model also predicts that the cavity should be filled 
with Tf. Ki 40, 000 K electrons giving rise to soft X-ray 
emissions. iLeahv et al.l (Il992[ l detected X-ray emitting 


plasma with Tg ~ 6 x 10® K towards the interior of the 
Gum Nebula, but we note that this coul d also be ex¬ 
plaine d by the wind-blown-bubble model of lWeaver et al. 


( 1977h . The wind-blown-bubble model also naturally ex¬ 
plains the ionised shell structure. 

The consensus in the literature to date favours the 
old SNR model of the ne bula, however, this i s not 
universally accepted (e.g., Ghoudhurv fc Bhatt 20091 : 
Urouhart et al. 2009l l. 


1.2. This work 

One way to differentiate between models of the Gum 
Nebula is to examine the density profile at the edge of the 
shell and the effect the nebula has on the magnetic field 
of the ISM. Supernovae and wind-blown-bubbles drive 
strong shocks into the ISM, compressing the gas at their 
leading edge. At the same time the gas inside the neb¬ 
ula may be ionised by the passing shock-front (in the 
supernova case) or by the central stars (in the case of 
wind-blown-bubbles) leading to a corresponding increase 
in electron density and magnetic field strength. Non- 
radiative shocks (for example in young supernovae less 
than ^ 20000 years old) expand adiabatically and we 
would expect to see a density compression factor V < 4 
at the edge of the shell. If the swept-up-shell has begun 
to cool radiatively (e.g., for snow-plough phase supernova 
remnants older than ^ 20000 years) then X can be much 
greater - up to several hundreds. Alternatively, if the 
bubble is due to a slow ionisation front moving into the 
medium, we would expect little compression and would 
measure X pe 1. 










































































































































The expansion of the bubble into the ISM should also 
imprint a clear signature on the Galactic magnetic field. 
The total field can be visualised as a superposition of an 
ordered large-scale component and a random small-scale 
component. The field lines are frozen into the gas, hence 
compression at the bubble edge can lead to an amplifica¬ 
tion of the field parallel the shock front. Faraday rotation 
is an especially sensitive probe of the field strength along 
the line of sight and this amplification is best observed 
as a rotation measure enhancement towards the limb of 
the shell. RMs also constitute an excellent probe of tur¬ 
bulence in the ISM. Unresolved random motions in the 
ionised gas can produce fluctuations in the random field 
that increase the scatter between adjacent RM samples 
and depolarise diffuse background polarised emission. 

In this work we combine point source measurements of 
rotation measures from background radio-galaxies, emis¬ 
sion measures (EMs) from Ha images and polarised 
2.3 GHz radio-continuum data to build a self-consistent 
picture of the Gum Nebula. Using a simple geometric 
model we derive the ambient electron density and mag¬ 
netic field strength. We fit for the compression factor 
in the shell, probe the geometry of the ordered Galactic 
field and shed light on the likely origin of the nebula. 


2. Datasets and Images 


We draw on data from seve ral publicly a v ailabl e sky 
surveys. We make use of the Tavlor et al. ( 2009ll RM 
catalogue, which is derive d from the 1.4 GHz NRAO 
VLA Sky Survey (NVSS, iGondon et al. 1996 1. We 
estimate emission me asures using the Southern H a 
Sky Survey (SHASSA, iGaustad et al. 2001 ; Finkbeiner 
2003h . dispersion measures (DMs) from the Australia 
National Telescope Fa cility Pulsar Catalogu^il (APG, 
Manchester et al.ll2005 1 and we examine the polarisation 
properties of the 2.3 GHz radio-continuum maps from the 
S-band Parkes All Sk y Survey (S-PASS, Garrettil 2011 


Garretti et al. 2013lJl . Below we introduce each of the 
surveys and describe the processing necessary to isolate 
the Gum Nebula from contaminating data. 

2.1. H a emission 

The n = 3 — 2 Balmer series H a recombination transi¬ 
tion of neutral atomic hydrogen is commonly used to 
derive EMs of ionised gas in the interstellar medium. 
EM is directly related to the electron density via 
EM = C rig dl, meaning that the intensity of Ha emis¬ 
sion can be used to estimate the line-of-sight electron 
density (see N3.3l for a full explanation). The So uthern H- 
Alpha Sky Survey Atlas ( Gaustad et al.ll20^ currently 
provides the highest spatial resolution (0fwhm = 6') 
coverage of the whole Gum Nebula in the Ha emission 
line. We use the r eprocessed SHASSA data published by 
Finkbeiner ( 20031) . who subtracted point-source emission 
from stars, corrected for imaging artifacts and calibrated 


the amplitude scale to the stable ze ro-point of the Wis- 
consin H-Alpha Mapper (WHAM, iHaffner et al.l 20031) 
survey. The H a image of the Gum Nebula is presented 
in Fig. [H The nebula describes a roughly circular shell of 
emission centred on the Galactic Plane. Below the mid¬ 
plane the structure of the H a data is very complicated, 
displaying arcs and filaments associated with overlapping 
Hii regions, SNR and other shells. The lower border 
of the Gum Nebula appears tenuous. Above latitudes 
b > 5° the structure of the Gum Nebula is much less 
confused. T he only obvious contamin ating feature is the 
Antlia SNR (iMcCullough et ah 2002), which overlaps at 
(/, b) w (268°, -1-12°). The northern arc of the Gum Neb¬ 
ula is particularly prominent, showing a sharp edge and 
a shell-like structure of width ^ 2°. 


2.1.1. Extinction correction 

H a emission is affected by extinction due to interven¬ 
ing dust along the line of sight, characterised by the op¬ 
tical depth T. If all of the dust responsible for the extinc¬ 
tion is in the foreground, then the intrinsic intensity /h„ 
is reduced by a factor e'^ to give the observed intensity 
dH„,obs- Since the location of the dust is unknown, the 
value /h„ = /hc ,obs can be considered a lower limit 
on the intensity (i.e., the maximum correction possible) 
and that of /h„ = /H„,obs an upper limit. If the dust 
is uniformly mixed with the source then the in trinsic in¬ 
tensity is given by /h„ = /H„,obs t / {1 — e~^) I Reynold^ 
h976all . 

In practice t may be determined from the extinction 
observed in the optical band, as it is related to the Eb-v 
colour by r = 2.44 x Eb-v ( Finkbeiner 2003 ). We 
corrected the H a data for ext i nction using the Eb-v 
map created by ISchlegel et al.l ( 1998h from the COBE 
(Cosmic Background Explorer) and IRAS (Infrared As¬ 
tronomical Satellite) surveys. These Eb-v maps pro¬ 
vide an estimate of the total column of dust in the 
Galaxy along the line-of-sight, however, towards higher 
latitudes it is reasonable to assume that most of the dust 
is nearby. Dust in the plan e has a scale-height of ^ 130 pc 
( Drimmel fc SDer^l200ll ) and at a latitude of b = 10°, 
sight-lines exit the dusty disk at a distance of ~ 800 pc. 

Upon inspection, images of /h„ produced assuming all 
dust is in front of the Gum Nebula appear over-corrected 
for prominent dust features. For example, a filament 
in the Eb-v map at (^, 6) = (250.4°, 14.5°) and a circu¬ 
lar feature at (Z, b) = (257.0°, 11.9°) turn from absorption 
features in /H„,obs to emission features in /h„. Thus, our 
best estimate for /h„ assumes that the dust is uniformly 
mixed with the H a emitting gas. At latitudes of 6 > 5° 
the values of r range over 0.20 < r < 0.93, co rrespond¬ 
ing to c orrections of 1.1 < r/(l — e~'^) < 1.5. iRevnold^ 
( 1976al) found r = 0.15 towards ( Puppis, corresponding 
to Ay = 0.19. This lower value of optical depth is consis¬ 
tent if we consider that the star lies just inside the front 
face of the Gum nebula. 


^http://www.atnf.csiro.au/research/pulsar/psrcat, VI.49, August- 
2014 










































































2.1.2. Galactic background 

Large-area H a images show that emission from dis¬ 
crete Galactic objects (e.g., Hii regions and supernova 
remnants) is superimposed on a diffuse background that 
rises to a peak at the Galactic mid-plane. As seen in 
Fig. [1] this is a particular problem for the Gum Neb¬ 
ula due to its large angular size and position straddling 
the plane. We have isolated the Ha emission from the 
nebula by estimating and subtracting a diffuse emission 
profile as a function of Galactic latitude. To calculate 
the profile we identified and masked-out all foreground 
objects in the image, took the minimum of the pixels 
in the longitudinal direction and smoothed the resultant 
profile to a resolution of ~ 1.6°. We initially created a 
background-corrected Ha map by subtracting a scaled 
version of this profile from each column of pixels in the 
original image. This simple scheme assumes that the 
background emission is constant with longitude across 
the 36° nebula; clearly not the case since Ha emission 
in the right hemisphere of the nebula is over-subtracted 
using this method. To further correct the background 
gradient, we fit an additional polynomial surface of or¬ 
der 2 to the residual large-scale emission. After back¬ 
ground correction, the brightness of the Gum Nebula’s 
shell varies between 30 R and 170 R away from the mid¬ 
plane (1 Rayleigh =10®/47rphotons s“^ cm“^ sr“^), com¬ 
pared to a mean background level of 2-7R. These val¬ 
ues are comparable with previous estimates made using 
pointed spectral observations capable of separating the 
Galactic and Gum Nebula components in velocity space 
( Reynold^ 1976bl l. The background-corrected Ha data 
are used in i i3.3l to estimate along the line of sight. 


2.1.3. Uncertainties 

The formal uncertainty in the Ha intensity is a 
quadrature sum of the intrinsic measurement uncertainty 
~ 0.3 R, and uncertainty due to the extinction 
correction Udust- As we do not explicitly know where 
the dust lies along the line of sight, we assume the worst 
case scenario and set the error to the likely range of 
correction values, typically Udust ~ 0.25. For values 
of iHct observed towards the northern arc of the Gum 
Nebula the absolute uncertainty is order 10 percent, or 
cr{lHa) ~ 12R in the shell. 

2.2. Rotation Measures 

Rotation measures of polarised background radio¬ 
galaxies provide a convenient method of measuring the 
line-of-sight magnetic field i?|| in local Galactic struc¬ 
tures, if the electron density ne and the distribution of 
the ionised gas are known. Each extragalactic point 
source is effectively at infinity and the measured differ¬ 
ence between RMs of adjacent sources is dominated by lo¬ 
cal changes in H|| or rie along the line-of-sight path. The 
effective RM contribution of a Galactic Hii region, for 
example, can be found by measuring the RMs of radio¬ 
galaxies behind the Hii region and subtracting an average 
off-source RM, determined from the radio-galaxies in the 


surrounding sky (e.g.. iHarvev-Smith et al.ll201lli . 

At the present time the best sampled and most ac¬ 
curate large-scale RM-grid covering a larg e part of the 


Gum Nebula is the cata logue created bv iTavlor et al. 


( 2009li from the NVSS ( Gondon et al.l Il996 b NVSS 
radio-continuum observations were conducted on the 
Very Large Array (VLA) at a frequency of 1.4 GHz and 
extend south to a declination of —40°. The original sur¬ 
vey combined simultaneous snapshot observations in two 
42 MHz-wide bands (1364.9 MHz and 1435.1MHz) into a 
multi-freque ncy synthesis i mage of the sky in Stokes /, 
Q, U and V. Tavlor et al.l ( 200^ reprocessed the NVSS 
visibility data into individual images of the bands and 
calculated two-channel rotation measures for 37,543 po¬ 
larised sources. 

2.2.1. RMs through the Gum Nebula 

Fig. [2]- left presents RMs from the Tavlor et al. ( 2009[ ) 
catalogue over-plotted on the H a image of the Gum Neb¬ 
ula. Although the lower left of the nebula is not cov¬ 
ered by the NVSS, there are RM measurements towards 
the upper right wi th source densities varying between 
1.2 and 6.6deg“^ ( Stil fc Tavlor 2007 1. It is easier to 
visualise RM-features iii the smoothed ma p created by 
Oppermann et all ( 2012h mainly from the iTavlor et al 


( 2009[l catalogue and plotted in Fig. [5]- right. Using the 


Ha and radio continuum maps as a guide, we identify 
structures in the RM map likely associated with Galac¬ 
tic objects. The RM-signature of the Gum Nebula is 
clearly different from the background and matches the 
morphology of the excited hydrogen gas well. The dis¬ 
tinctive upper arc of the nebula displays consistently pos¬ 
itive RMs, which decrease in magnitude towards the ge¬ 
ometric centre, reminiscent of a limb-brightened shell. 
The net positive RM signal in the upper Gum Nebula 
implies a coherent magnetic field on scales of ~ 260 pc, 
the projected diameter of the nebula at the adopted dis¬ 
tance oi d = 450 pc. The region along the Galactic mid¬ 
plane (l&l < 5°) is highly confused, containing several 
Hii regions and supernova remnants (see Fig. [T|), while 
the lower part of the nebula also contains the IRAS Vela 
Shell, which may be a separate foreground object. In 
contrast, the upper part of the Gum Nebula appears rel¬ 
atively free of obscuring objects and we focus on this 
region for the remainder of the paper. The grey wedge- 
shaped box in Fig. [2] outlines the data selected for anal¬ 
ysis. The selected region includes the upper arc of the 
nebula, the interior above 6 = 5° and a section of off- 
source RMs outside of the nebula’s border. 

2.2.2. Isolating the RM-signature of the Gum Nebula 

Any analysis of the RM data relies on isolating the 
RM-signature of the Gum Nebula from discrete regions of 
magneto-ionic material along the line of sight (e.g., over¬ 
lapping supernova remnants or H ii regions) and from 
the bulk of the Galaxy in the background. We have 
identified three other sets of discrete Galactic objects 
towards the upper G um Nebula that hav e Faraday rota¬ 
tion signatures in the Tavlor et al. ( 20091) RM-catalogue. 
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Fig. 2.— Left: The RM-catalogue of iTavlor et alJ (|2009l ) plotted over the Ha map of the Gum Nebula (|Finkbeinerll200 . Red 
circles indicate positive RMs, while blue indicate negative and their diameter is proportional to |RM|. The solid Ha contour at 
a level of 25 R def i nes th e outline of the nebula. Right : Prominent RM-features are easier to visualise in the map produced by 
lOppermann et ^ l|2012l ^ using the lTavlor et alJ (|2Q09l ) catalogue. Polygons and lines annotate significant features. We restrict 
our analysis to the upper arc region, inside the solid grey line. 
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Fig. 3.— Left: Positions of polarised extragalactic sources selected for analysis from the lXavlor et al.! (l2009h catalogue. Red 
and blue circles indicate positive and negative RMs, respectively, and the radius of each circle is scaled to the absolute value of 
the RM. The green dashed lines enclose regions where sources have been excised from the sample because they likely probe a 
contaminating region, or are obvious outliers (> 3 (t). Right: RMs plotted as function of angular offset from the centre of the 
nebula. The formal uncertainties in the catalogue are represented by vertical grey lines (scaled x 2 for clarity) and each point is 
colour-coded for Galactic Longitude. The solid black histogram represents a median-binned version of the plot, with bin-sizes 
of 0.6°. The solid grey histogram illustrates the level of the RMs exterior to the nebula before the discrete object correction 
was applied. Vertical black error-bars illustrate the la scatter in each bin. The magenta line in both panels indicates the outer 
boundary of the nebula. In both panels the RMs are plotted assuming a flat background, i.e., no correction has been made 
for a large-scale RM-gradient due to the Galaxy in the background (or foreground). The effect of subtracting different model 
backgrounds is illustrated in Fig. |4]and discussed in 

















































The Antlia SNR (jMcCullough et alJ 120021 annotated in 
Fig. [2]) is adjacent to the nebula on the upper left. The 
edge of the SNR is traced by sources with RM values 
10 — 15radm“^ more positive than their surroundings. 
This excess RM signal is comparable to the formal er¬ 
ror in the catalogue 12radm“^) and is negligible 
compared to RMs through the rim of the Gum Neb¬ 
ula 300radm“^). RMs through the interior of the 
Antlia SNR are consistent with the background, except 
for a patch directly bordering the Gum Nebula, which 
is more negative than the large scale background by ap¬ 
proximately —30.0radm“^. This patch lies inside our 
selection box so we subtract this offset from RMs inside 
the patch to correct the catalogue. The second obvi- 
ous feature in t he da ta is a pair of shells discovered by 
lacobelli et al. ( 2014l l in 2.3 GHz radio continuum data, 
seen in polarisation and lying to the upper-right of the 
Gum Nebula. The border of the shells seen in the radio 
data corresponds exactly to the morphology of a neg¬ 
ative patch of RMs in Fig. [2l We again correct the 
catalogue by subtracting the median background offset 
(—77.2radm“^) from RMs inside the shell boundaries. 
The final feature of note is a ‘stalk’ of positive RMs ex¬ 
tending from the centre of the upper arc to higher Galac¬ 
tic latitudes. Th e ‘stalk’ has a coun t erpar t in Hi emis¬ 
sion identified bv iRevnoso fc Dubnerl (1993), lies above a 
hole in the H a image and is hypothesised by the authors 
to be a ‘blowout’ in the shell wall leading to ionised gas 
streaming into the Galactic halo. Modelling and sub¬ 
tracting the signature of this feature is beyond the scope 
of this work, so we simply mask off the RM data within 
its boundary. Fig. [3l- left presents the corrected RM cat¬ 
alogue within the selection box, plotted over the Ha 
image. The azimuthally-averaged RM-profile is shown 
in Fig. right). The black histogram shows a version 
of the profile binned in 0.6° increments. Outside of the 
nebula border (offsets > 22°) the RMs are relatively con¬ 
stant but rise rapidly to a peak just inside the border. 
At smaller offsets the RM values fall slowly, approaching 
an interior level that is higher than the background. The 
grey histogram, shows the same binned profile prior to 
correcting for discrete contaminating sources. Note that 
the RM data shown here have not yet been corrected for 
large-scale gradients due to diffuse thermal electrons dis¬ 
tributed throughout the bulk of the Galaxy in the back¬ 
ground. 

The distribution of electrons within the Galaxy and the 
strength, and geometry, of the ordered Galactic magnetic 
field result in a unique pattern of RMs over the whole sky. 


In the all-sky RM map compiled by lOopermann et al 


(|2012ri the dominant signal is quadrupolar in shape, with 
negative RMs above and positive below the Galactic 
mid-plane in the vicinity of the Gum Nebula. In re¬ 
cent years several authors have modelled the ordered 
Galactic field by combining data fro m extra-galactic 


RMs and radio-synchrotron emission (ISim^^_aI 


.lansson et al. 

2OO9I 

Sun & Reich 

20IC 

: Mao et al. 

2010; 

Jaffe et al.l 

201 nt Van Fck et al .1120111: 

.lansson & Farrar 


2008 


toroidal field in the halo. Within the Galactic disk 
the magnetic field and thermal electron density follow 
the spiral arms, increasing towards the mid-plane, lead- 
ing to steep gradients in RM at low Galactic latitudes 
jSimaixFNormandm Kr onberd 1980l : iGordes fc Lazio 
200^ Gaens]er^t_^ 2008ll . The sparse sampling of the 


Tavlor et al.l (|2009ll RM-catalogue and confusion towards 


the mid-plane mean that accurately removing the large- 
scale RM-signal due to the Galaxy is challenging. We 
initially attempted to fit a 2D polynomial surface to the 
off-source RMs, but this proved to be highly unreliable 
in practice. Instead we consider two classes of poten¬ 
tial RM backgrounds. In the first case we assume a 
simple flat background at the median level of the se¬ 
lected RMs outside the boundary of the Gum Nebula: 
—26.4radm“^. This assumption is the simplest correc¬ 
tion possible and consistent with the high-latitude RM- 
data. However, we know that the volume-averaged elec¬ 
tron density decays exponentially with height above the 
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F ig. 4.— T o p: Sa mple profiles at I = 259° from the models 
of ISun et al.l (120081 1 and ijansson fc Farrar! 1 I 2 OI 2 I I illustrating 
the predicted RM-gradient with Galactic latitude. Bottom: 
The light-grey histogram is the azimuthally averaged RM- 
profile of the Gum Nebula uncorrected for discrete contami¬ 
nating sources or large-scale RM background (also plotted in 
Fig. El . The black histogram is the equivalent with discrete 
corrections applied and a flat background of —26.4radm“^ 
subtracted. The red (dashed) and green (dot-d ashed) his¬ 
togram s ill ustrate the effect of subt racting the 2-D ISun et ahl 
(I 2 OO 8 II and ijansson fc Farrarl 1 I 2 OI 2 I) RM- models, respe c tively , 
from each point source select from the I Tavlor et al.1 (l2009l l 
catalogue. Note that the discrete corrections have previously 
been applied. In both cases the large-scale gradient correc¬ 
tion acts to decrease the RMs in the interior of the nebula 
relative to the exterior. 


































































































mid-plane ( Cordes fc Lazio 2002 : l^ensler et alJ 2008 ). 


thus t he R Ms must decrease corresp ondingly. iSun et al 


( 2008h and Jansson & Farrar ( 2012h have modelled the 
large scale RMs distribution of the Galaxy starting from 
the NE2001 electron density dis tribution and applyin g 


the scale height corrections of Gaensler et al. ( 2008ll . 


Both models have similar latitude profiles, illustrated for 
I = 258° in Fig.|4]- top. The bottom panel of Fig.|4]shows 
the effect of subtracting each 2-D model from the selected 
RM data-points. The resulting azimuthal profile is essen¬ 
tially the same for both models, but offset in RM as the 
models differ in their absolute calibration. Both models 
act to decrease the value of RMs towards the interior of 
the nebula. Neither model correctly predicts the abso¬ 
lute zero-level exterior to the Gum Nebula, likely because 
the best-fitting models are constrained over the whole sky 
and by other, sometimes contradictory, datasets. The au¬ 
thors also had limited knowledge of lo cal contam i nating 
obje cts. We adopt the RM-m odels of LSun et al. (2008) 
and I.Tansson fc Farrail ( 2012[l as the best available es¬ 
timates of the large-scale background variation in RM 
and apply offsets of —64.0radm“^ and —36.6radm“^, 
respectively, so as to correct their calibration to the zero- 
point exterior to the Gum Nebula. In our analysis of the 
RMs through the Gum Nebula we compare the results 
derived with each of these backgrounds separately. 

2.3. 2.3 GHz radio continuum 

Radio continuum at centimetre wavelengths traces gas 
emitting via both synchrotron and thermal processes. If 
information on the polarisation state of the radiation is 
available, analysis of the Stokes Q and U parameters can 
constrain conditions in the gas along the line-of-sight, 
e.g., depolarisation due to a fluctuating component of 
the magnetic field. 

The S-band Polarisation All Sky Survey (S-PASS) has 
imaged the entire southern sky (Dec. < —1°) in polari¬ 
sation at a frequency of 2.3 GHz. The observations have 
been conducted with the Parkes Radio Telescope, NSW 
Australia, a 64-m telescope operated by GSIRO Astron¬ 
omy and Space Science. A descr iption of S-PASS ob ¬ 
serv ations and ana l ysis is given in Garretti et al. ( 2010ll 
and Garretti et al.l (2013b). Here we report a summary 
of the main details. The standard S-band receiver of 
the observatory (Galileo) was used with a system tem¬ 
perature Tsys = 20 K, beam width FHWM= 8.9' at 
2300 MHz and a circular polarisation front-end ideal 
for linear polarisation measurements with a single-dish 
telescope. Data have been detected with the Digital 
Filter Banks mark 3 (DFB3) with full Stokes capabili¬ 
ties recording the two autocorrelation (RR and LL) and 
the complex cross-correlation products of the two circu¬ 
lar polarisations (RR, LL, LR, RL*). Flux calibration 
was done with PKS B1934—638, secondary calibration 
with PKSB0407—658 and polarisation calibration with 
PKS B0043-424. Data were binned in 8 MHz channels 
and, after RFI flagging, 23 sub-bands were used, cov¬ 
ering the ranges 2176-2216 and 2256-2400 MHz, for an 
effective central frequency of 2307 MHz and bandwidth 


of 184 MHz. 

The observing strategy is based on long azimuth scans 
taken towards the East and the West at the elevation of 
the south celestial pole at Parkes (EL= 33°) to realise 
absolute polarisation calibration of the data. Final maps 
are convolved to a beam of FWHM= 10.75'. Stokes I, 
Q, and U sensitivity is better than 1.0mJybeam“^ per 
beam-sized pixel everywhere in the covered area. De¬ 
tails of scanning strategy, map-making, and final maps 
ob tained by binning all frequency channels are presented 


Garretti et al. ( 2010ll and Garretti et al. (2015, in 


preparation). The co nfusion limit is 6 mJy in Stokes I 
( Garretti et al.l [2013all and much lower in polarisation 
(average polari sation fraction in c ompact sources is lower 
than 2 percent. iTucci et al. 2004h . The in strumental po¬ 


larisa tion leakage is 0.4 percent on-axis (jGarretti et al 


2010[j and less than 1.5 percent off-axis. For diffuse emis¬ 
sion, the latter is generally not important because of 
cancellation effects at scales larger tha n the beam (e.g., 
Garretti et al. 2004 lO’Dea et al. 2007 1. 


2.3.1. The Polarised Signature of the Gum Nebula 

Fig. [5] presents the 2.3 GHz radio-continuum image of 
the Gum Nebula in Stokes /, Q, U and polarised inten¬ 
sity P. The morphology of the nebula in total intensity is 
broadly similar to the Ha map presented in Fig. [TJ im¬ 
plying that the radio- and optical-emission are coming 
from the same gas. When viewed in P, Q and U, the up¬ 
per shell of the nebula is seen to depolarise background 
emission in a ~ 2° wide arc. This band of depolarisation 
is set against the smooth Galactic background, visible in 
the upper-right quadrant of the image, above b = 12°. At 
high latitud es the background is also depolarised by two 


thin shells (jiacobelli et al.l 12014 see 112.2.11 and Fig. [2]) 


and by the Antlia supernova remnant in the upper-left 
quadrant. The shell of the Antlia SNR is similarly char¬ 
acterised by a band of depolarisation that overlaps the 
Gum Nebula at (l,b) ps (268°,-|-13°). The interiors of 
the Gum Nebula and Antlia SNR appear fractured, ex¬ 
hibiting patches of homogeneous polarised intensity in¬ 
terspersed with depolarised ‘canals’. The Vela supernova 
remnant at (I, b) = (267°, —3°) is the brightest object in 
the field, (I & P), while the rest of the Galactic plane 
is seen as a mix of polarised foreground and depolarised 
background emission. 

2.4. Pulsars 


Th e ATNF Pulsar Gatalogue (APG. [Manchester et al 


2005h collates the properties of more than 2300 rotation- 
powered pulsars and is continually revised as new dis¬ 
coveries are made. Of primary interest to us are the 
dispersion measures (DMs), RMs and distances to the 
pulsars. These parameters can be combined with EMs 
and a geometric model to derive the average electron den¬ 
sity, filling factor and magnetic field strength along the 
line-of-sight. 

We utilise version 1.49 of the APS, which lists 158 pul¬ 
sars within a 30 ° radius of the kinematic centre of the 


Gum Nebula (I = 261.0°, b = —2.5°: IWoermann et al.l 
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Fig. 5.— 2.3 GHz radio-continuum maps of the Gum Nebula from the S-PASS project. The data are calibrated in Janskys and 
may be converted to a main-beam brightness temperature scale in Kelvin by multiplying by 0.55. The border of the nebula 
traced by the green line is the same as in Fig. [5] 


200lh . Of these we have chosen 35 for analysis, most 
of which lie within the upper Gum region. Fig. [ 5 ] and 
Table [T] present this sample, which contains all pulsars 
above b — 2° and a handful below, chosen because they 
have accurately determined distances or lie on uncon¬ 
fused sight-lines adjacent to the Gum Nebula. 

Accurate distances to pulsars are difficult to obtain: 
a handful of precise values have been calculated via an¬ 
nual parallaxes and these are limited to relatively nearby 
pulsars (< 3kpc). Kinematic distances accurate to 

1 kpc can be derived for some pulsars associated with 
Hi absorption, while the distance to pulsars located in 
globular clusters can be estimated to 15% reliabil¬ 
ity via analysis of colour-magnitude diagrams. Most 
of the pulsars detected towards the Gum Nebula de¬ 
fault to a distance derived from the dispersion measure. 
Such DM-distances are often highly inaccurate because 
they rely on a model of the Galactic free-electron dis - 
tribution ( Cordes fc Lazio 20021 : iTavlor fc Corded 1993 1. 


which was itself created in part using the iTavlor et al 
I I 993 II pulsar catalogue. Distances are particularly ill- 
determined t owards the Gum Neb ula, which was in¬ 
cluded in the ICordes fc Lazio (2002) model as a pair of 
overlapping spheres of diameter 50 pc. I t is not clear that 
this is an improvement on the older ITavlor fc Gordes 
1 1 993(1 model which treated the nebula a simple Gaus¬ 
sian of full width half m aximum ('FWHM) 50 p c trun¬ 
cated at r = 130 pc. The Cordes fc Laziol (2002) model 
does, however, account for the scatter broaden i ng T^r 
which was measured bv iMitra &: Ramachandr^ ( 2001 1 
for 40 pulsars between 250° < I < 290°. They found 
that Tsc was greater than expected for a smooth Gaus¬ 
sian, implying a more inhomogeneous distribution of rie- 
Based on the observed scattering they concluded that 
pulsars in the vicinity of the G um Nebula should be 2 - 3 
times closer than predicted bv ITavlor fc Corded ( 1993[l . 
Within the area of the Gum Nebula only four pulsars 
have both DM measurements and accurate distances. 
































































Table 1: Properties of selected pulsars towards the Gum Nebula. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Name 

1 

b 

DM 

RM 

D 

Notes 


(deg.) 

(deg.) 

(cm”^ pc) 

(radm“^) 

(kpc) 



Pulsars towards the upper Gum Nebula: 


J0758-1528 

234.464 

07.224 

63.327 ±0.003 

55 ±7 

3.72 

Out 

J0818-3049 

249.983 

02.908 

133.7±0.2 

- 

4.17 

Gum 

J0820-1350 

235.890 

12.595 

40.938 ±0.003 

-1.2±0.4 

1.9+ 

Out 

J0828-3417 

253.965 

02.561 

52.2±0.6 

59 ±3 

0.5 

Gum 

J0838-2621 

248.807 

08.981 

116.9±0.1 

86± 13 

4.6 

Gum 

J0846-3533 

257.190 

04.710 

94.16±0.11 

144 ±8 

1.4 

Gum 

J0855-3331 

256.847 

07.517 

86.635 ±0.016 

165 ±10 

1.2 

Gum 

J0900-3144 

256.162 

09.486 

75.702 ±0.010 

- 

0.8 

Gum 

J0904-4246 

265.075 

02.859 

145.8 ±0.5 

284 ± 15 

4.4 

Gum 

J0908-1739 

246.119 

19.850 

15.888 ±0.003 

-31±4 

0.6 

Out 

J0912-3851 

263.165 

06.584 

70 ± 1 

- 

0.6 

Gum 

J0923-31 

259.697 

13.003 

72 ±20 

- 

1.0 

Gum 

J0932-3217 

261.277 

14.069 

102.1 ±0.8 

- 

3.8 

Gum 

J0934-4154 

268.361 

07.411 

113.79 ±0.16 

- 

3.2 

Gum 

J0941-39 

267.795 

09.904 

78.2±2.7 

- 

1.3 

Gum 

J0945-4833 

274.199 

03.674 

98.1 ±0.3 

- 

2.7 

Out 

J0952-3839 

268.702 

12.033 

167 ±3 

- 

8.4 

Gum, Ant 

J0959-4809 

275.742 

05.418 

92.7± 1.2 

50 ±6 

3.0 

Out 

JlOOO-5149 

278.107 

02.603 

72.8 ±0.3 

46 ±9 

2.3 

Out 

J1003-4747 

276.037 

06.117 

98.1 ±1.2 

18 ±4 

3.4 

Out 

J1012-2337 

262.131 

26.377 

22.51 ±0.09 

52 ±9 

1.3 

Out 

J1032-5206 

282.354 

05.128 

139 ±4 

- 

4.3 

Out 

J1034-3224 

272.050 

22.117 

50.75 ±0.08 

-8± 1 

4.7 

Ant 

J1036-4926 

281.518 

07.727 

136.529 ±0.010 

-11±6 

8.7 

Out 

J1045-4509 

280.851 

12.254 

58.166 ±0.001 

92 ±1 

0.23+ 

Ant 

J1057-4754 

284.007 

10.739 

60 ±8 

- 

3.0 

Ant 

J1105-43 

283.511 

14.886 

38.000 ±0.001 

- 

2.2 

Ant 

Pulsars below b — 2° with accurate distances: 




JlOOl-5507 

280.226 

00.085 

130.32 ±0.17 

297 ± 18 

0.3+ 

Out 

J0737-3039B 

245.236 

-4.505 

48.920 ±0.005 

112.3 ± 1.5 

1.1 + 

Out 

J0738-4042 

254.194 

-9.192 

160.8 ±0.7 

12.1 ±0.6 

1.6+ 

Gum 

J0742-2822 

243.773 

-2.444 

73.782 ±0.002 

149.9±0.1 

2.0+ 

Out 

J0835-4510 

263.552 

-2.787 

67.99 ±0.01 

31.4±0.1 

0.28+ 

Gum 

J0837-4135 

260.904 

-0.336 

147.29 ±0.07 

135.8 ±0.3 

1.5^ 

Gum 

J0908-4913 

270.266 

-1.019 

180.37 ±0.04 

10.0± 1.6 

1.0+ 

Gum 

J0942-5552 

278.571 

-2.230 

180.2 ±0.5 

-61.9 ±0.2 

0.3+ 

Out 


Note.— This table presents the properties of the pulsars displayed in Figure^ Pulsars marked with a f in column (6) have independent 
distance measurements; other entries default to the DM-derived distance. The code in column (7) notes whether the pulsar falls on a 
sightline towards the Gum Nebula (Gum), towards the Antlia SNR (Ant) or outside of the border of either object (Out). 


The Vela pulsar (J0 835—4510) is known to be at a dis¬ 
tance of 287^17 pc I Dodson et al. 20031) . placing it just 
inside the front wall of the nebula. The remaining three 
(J0738-4042, J0837-4135 and J0908-4913) lie at dis¬ 
tances greater than Ikpc (see the bottom of Table [T]), 
behind the Gum Nebula. 


3. Analysis 

Our analysis aims to answer the questions: What is 
the likely origin of the Gum Nebula? and What are the 
magnetic properties of the nebula and how do they af¬ 
fect ambient conditions in this part of the Galaxy? To 
address these questions we construct a simple model of 
the nebula as an ionised shell situated in the near-held. 
We present the model below and explain the maximum- 
likelihood method used to ht the model to RMs on the 
sky. The model assumes a uniform density distribution 
plus a jump in ionisation fraction from 0 to 100 percent 
within the shell of the Gum Nebula, which we derive 
from the H a data and include as a prior in our Htting 
procedure. The resulting hts will be presented in 2) 


3.1. Rotation Measures as Magnetic Probes 

Faraday rotation causes the polarisation angle of a 
linearly polarised wave traversing a magnetised ionised 
medium to rotate by an angle Nip. The change in polar¬ 
isation angle is given by 

Ai/; = RMA^ rad, (1) 

where RM is the rotation measure in radians m“^. The 
observed RM depends on the line-of-sight component of 
the magnetic held R|| (in ^G), the thermal electron den¬ 
sity Ue (in cm“^) and the path length dl (in pc) according 
to 

nobs 

RM = 0.81 / UeB^^dl radm“^. (2) 

J src 

Note that the integral in Equation [2] is taken from the 
source of the polarised emission to the observer, so that 
a positive RM indicates an average magnetic held point¬ 
ing towards the observer. If the ionised material along 
the line-of-sight contains clumps of uniform Ue threaded 
by the same i3||, then the medium is characterised by a 
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Fig. 6.— Selection of pulsars with known DM values towards the Gum Nebula. All pulsars with 6 > 2° are shown alongside a 
selection of pulsars with accurate distances below b = 2°. Red dot-dashed lines outline confusing parts of the Galactic plane and 
pulsars in these regions have been omitted, ‘o’ symbols represent pulsars with both DM and RM measurements, and ‘-I-’ symbols 
pulsars with only DM measurements. Each pulsar is annotated with its DM in pccm“® (top-left number, coloured green), RM in 
radm“^ (bottom-left number, coloured red or blue signifying positive or negative RM, respectively), independently determined 
distance in kpc, where known (top-right number, coloured black) and EM in pccm~® (bottom-right number, coloured magenta). 


volume filling factor /. Equation [5] becomes 

RM = 0.81neB||/L radm-^, (3) 

where L is the total path length through the ionised 
medium and f L is known as the occupation length. L 
can generally be estimated from the geometry of the ob¬ 
ject under consideration (e.g., a slab, sphere or shell). 


3.2. A Near-Field Magnetic Bubble Model 

Models of RMs through spherical ionised shells have re¬ 
cently been used to derive magnetic properties of Galac¬ 
tic SNRs and Hii regions, and to probe the magnitude 
and orientation of t he or de red Galactic magnet i c field , 
e.g. Savage et afl (12013 ). Harvev-Smith et ^ ( 2011 1. 
Whiting et al.l ( 200^ and Kothes &: Brown ( 2009 1. 
These phenomena ionise their surroundings and illu¬ 
minate the ambient magnetic field via Faraday rotation. 
As they expand into the ISM they may also compress 
the field, imprinting a specific signature on the rotation 
measures. Previous investigations have focused on dis¬ 
tant objects (> Ikpc) whose small angular diameters 


(< 5°) mean that they intercept fewer RM sight-lines 
compared to the nearby Gum Nebula. Because these 
bubbles lie in the far-held, their RM proHles may be 
integrated in azimuth under the assumption of spheri¬ 
cal symmetry. However, for a near-held bubble like the 
Gum Nebula, the sign and shape of the RM-prohle can 
vary with azimuth, depending on the orientation of the 

ordered magnetic held. _ 

In a similar way to lWhiting et al. ( 2009ll and Savage et al 
(|2013l) , we model the Gum Nebula as a spherical ionised 
shell of radius R and thickness dr, threaded by a uni¬ 
form, parallel magnetic held Bq. We assume that the 
electron density Ue is constant within the shell and zero 
elsewhere (i.e., the background has been removed as in 
il2.1l and the electron density in the interior of the shell 
is negligible). The observer is located in the near-held 
and the magnetic held lines make an angle 0 to the 
plane of the sl^ in the direction of the bubble centre. 
The tilt anglqj of the magnetic held is assumed hxed 


^The tilt angle is the angle the magnetic field makes to the Galactic 
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Fig. 7.— Grid of models showing how changes in the angle of the magnetic field 0 and the compression factor X affect the 
distribution of RMs across a simple ionised shell. The geometry of the model is described in the Appendix (see Fig. A 1.) and 
the parameters are set to: D = 450 pc, (jiouter = 22.7°, dr = 25.0 pc, rie = 1.7 cm Bq = 8.6 and / = 0.5. Each model 
is presented in two panels: the upper panel presents a map of RM in offset Galactic coordinates and the lower panel displays 
radial RM profiles extracted over a range of angles T to the Galactic plane (0° < T < 180°). 

















































































































































along the y-axis, representing the Galactic plane, and 
the angle ( describes the orientation of the sight-line to 
the yz-plane. The edge of the shell subtends an angular 
radius of (/)outer = sm~^(E/D), where D is the distance 
from the observer to the geometric centre of the shell. A 
full description of the adopted geometry can be found in 
Appendix]^ including a detailed schematic. 

If the shell is expanding supersonically into the ISM 
then the gas, and hence the magnetic field, will be com¬ 
pressed at the external boundary. If the expansion has 
slowed to sub-sonic speeds then the gas will simply move 
out of the way. To model the compression (or lack of) we 
assume the electron density behind the expansion front 
is given by rie = X uq, where X is the compression factor 
and no is the electron density in the ambient medium. At 
each point on the sphere the component of the magnetic 
field tangent to the shell {B±) is amplified by X while 
the normal component (Bn) is unaffected. The contribu¬ 
tion to the observed magnetic field by one hemisphere is 
simply the vector sum of X B± and Bn projected along 
the line of sight. The measured RM is then proportional 
to the sum of the ingress (far hemisphere) and egress 
(near hemisphere) components. Equation |3] becomes a 
function of polar coordinate (<p,C)^ compression factor 
X, electron density Ue, magnetic field strength Bq and 
the angle of the magnetic field to the plane of the sky 0 


RM = 0.81/ 


/i3|l(</.,C,Ro.e,X) 

V hG 



^ L((f),dr) '^ 

(4) 


The model implicitly assumes that the same i3|| applies 
everywhere along the half-chord between the outer sur¬ 
face and mid-plane of the shell (but different for ingress 
and egress). This assumption will only be realistic for 
a thin shell; a more sophisticated analysis is outside the 
scope of this paper. The model also assumes a constant 
value for the electron density within the shell and, be¬ 
cause Ue, Bios and / are degenerate in Equation [31 the 
electron density and filling factor must be estimated from 
independent data, if possible (see H3.31 below). 

Fig. [7] presents a grid of near-held bubble models, il¬ 
lustrating how the RM observed on the sky changes as 
0 and X are varied. In the far-held case, the rotation- 
measure prohle is spherically symmetric, i.e., constant 
with (the angle between the prohle sampling line and 
the Galactic plane, see Fig. dH). However, in the near- 
held case the prohle shape depends on both (/ and 0 (the 
orientation of the ordered magnetic held vector). As 0 
increases from 0° (Bq in the plane of the sky) to 90° (Bq 
pointing towards observer) the large scale distribution 
of RMs changes from being anti-symmetric to symmet¬ 
ric around the central latitude axis. Increasing X leads 
to an overall increase in the magnitude of the RM and a 
large difference in RM between the centre and inner-edge 
of the shell. The shape of the prohle edge also becomes 
more rounded because of limb-brightening, although this 
would be much more noticeable in far-held bubbles. 

In tl3.41 below, we present a maximum-likelihood 
method used to ht the model to rotation measures from 


disk at the position of the nebula 


the Tavlor et al. ( 2009ll catalogue. We incorporate esti¬ 
mates of He from Ha data as a prior, assuming Gaussian 


errors. 


3.3. 


from H a data 


From Equationj^we see that rie and H|| are degenerate 
and cannot be determined individually from observations 
of RM alone. However, a separate observation of the 
emission measure provides an independent estimate of 
the electron density along the line of sight. EM is related 


to rip via 


(5) 


EM = / n/,dl pc cm °. 

Jo 

Assuming the same clumpy medium and geometry as 
presented in Section l3. 1 1 this becomes 


EM = rip f L pc cm 


-6 


( 6 ) 


with hlling factor / and path length L as before. EM may 
be calculated directly fr om the int e nsity of the H a (3 —2) 
line via the equation of iRevno Idi (1 198811 


EM = 2.75 


Tp 


104 K 


0.9 


Ihp.K 

R 


cm ®pc. 


(7) 


where Tp is the electron temperature in K, /h„ is the 
intensity of the Ha emission in Rayleighs and K = 
T / (I — e~'^) is a correction term to account for dust ex¬ 
tinction between the H a emission and the observer (see 
Section [2.1.111 . 

Most estimates of Tp for the Gum Nebula within the 
literature vary between 6500 K and 1150 0 K. Electron 
temperatures derived bv iRevnoldsI (1976a) from a com¬ 
parison of Ha to [Nil] linewidths are consistent with 
a uniform temp eratur e of 11 300 K throughout the neb¬ 
ula. However, Vidal! ( 1979ll determined the electron- 
temperature to be Tp = 6500 K from existing optical 
emission-line data. We adopt a uniform value of 8000 K 
for this analysis. 

Fig. [8]- top-left presents the emission m easure map 
of the upper Gum region derived from the iFinkbeineil 
1 2003h H a data using EquationjTl Azimuthally averaged 
EM values peak at 220pccm“® (see Fig. |5|- bottom-left), 
falling to ~ 80 pc cm“® in the interior and < 30 pc cm“® 
outside the nebula. Our value s are larg e ly con sistent 
with those of previous authors. IRevnoldsI (Il976all mea¬ 
sured the EM via pointed H a spectral observation and 
found it varied from 100pccm“® interior to the nebula 
to 240 pc cm~^ in th e upper arc. Scans across the region 
bv IRevnoldsI (1976a) at I = 240° determined the Galactic 
background to be < 10pccm“® rising to 28pccm“® at 
the mid-plane. 

Assuming the geometry L((j), dr) from the best-fit to 
our shell model (presented in [JH below) and best-fit 
filling factor / = 0.3 we use Equations [6] and IA2I to 
derive the electron density rip inside the clumpy shell. 
The electron density map and azimuthally averaged pro¬ 
file are presented in Fig. [8]- right. The fitting proce¬ 
dure takes the value of rip as a prior (see N3.4I) alter¬ 
ing the most likely shell geometry and filling factor /, 
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Fig. 8.— Left: EM map and azimuthally averaged profile of the upper Gum region. The individual grey points in the bottom 
panel show the EM values at the positions of the RM samples while the black points with error bars show the azimuthal average 
and standard-deviation in 0.2 degree bins. Right: Electron density map and azimuthal averaged profile derived assuming the 
best-fit parameters in Tabled The magenta line shows the edge of the shell model. 


requiring a new estimate of ng. To correct for this in¬ 
consistency we re-calculated Ue using the new L{(p,dr) 
and /, and iterated over the fitting loop until all values 
converged. The final electron density was determined 
to be Up, = 1.4 ± 0.4cm“^, which compares well with 
iRevnolda ( 1976tj) who found Ue ~ 1 cm“^ in the fainter 
parts of the nebula, assuming fixed physical parameters 
{4'outer = 18°, D = 450pc and 15 < dr < 30pc). 

3.4. Maximum likelihood analysis 

We use a Markov Chain Monte Carlo (MCMC) al¬ 
gorithm to fit the model shell to the RM data in the 
upper Gum region. The de-facto algorithm for per¬ 
forming MCMC fitting is the Metropolis-Hast ings al¬ 
gorithm ( Metropolis et al. 1953 : Hastingslll970ll . which 
randomly samples over parameter space, accepting or 
rejecting models based on their likelihood C (i.e., the 
probability of the data given the model parameters). 
New positions with greater C than previously are al¬ 
ways accepted, while those with smaller C are occasion¬ 
ally accepted. Ou r code makes use of the e fficient afhne- 
invariant sanipler ( Goodman fc Wei ^LoiOi) implemented 
in the emceeH pyt/ion module bv iForeman-Mackev et ahl 
(l2m.3h . EMCEE controls a number of parallel samplers, 
referred to as ‘walkers’, each of which corresponds to a 
vector of free parameters within the model. The walk¬ 
ers are initialised to a point in n-dimensional parame¬ 
ter space and are iteratively updated to map out the 


®http://dan.iel.fm/emcee/ 


probability distribution. At each iteration the likeli¬ 
hood is calculated assuming Gaussian errors according 
to C = e~^ where is the standard chi-squared 
goodness-of-ht statistic. If priors with measured uncer¬ 
tainties exist for any model parameter, we incorporate 
them into the likelihood calculation by summing their 
chi-squared values = xlyodei + E* Xpriori- 

The fitter is started by generating 300 walkers ini¬ 
tialised to random values of the free parameters. The 
MCMC code is initially run for 400 ‘burn-in’ iterations 
to allow the walkers to settle in a clump around the peak 
in likelihood space. The fitting routine is then run for 
10000 iterations to produce a well-sampled likelihood dis¬ 
tribution. We determine the best fitting model from the 
mean of the marginalised posterior distribution for each 
free parameter. The ilcr uncertainties are calculated 
as the fractional positions at 1 — erf(l(T) = 0.1572 and 
erf(lcr) = 0.8427 on the normalised cumulative distribu¬ 
tion. Our results are presented in M.2I 

3 . 4 . 1 . Scatter in RM as a hyperparameter 

The median measurement uncertainty on the selected 
RM data is (t(RM) = 12radm“^, considerably smaller 
than the scatter evident in Fig. [3]- right. The error-bars 
reflect only the uncertainty in the measurement and do 
not take into account systematic scatter, e.g., due to fluc¬ 
tuations in Rji or Ue on scales much smaller than the 
sampling grid, or systematic errors in RM determina¬ 
tion. We characterise this additional variation using a 
term (5(RM) added in quadrature to the measurement 


















































































































uncertainty. This new scatter term is included in the 
model as a free parameter, however, it is treated slightly 
d ifferently when calc ulatin g the likelihoo d funct ion. 


nerently wnen calc ulatin g tne likelinoo d tunct ion. 
Lahay et al. ( 20001 ) and iHobson et ahl ( 20021 ) present 


a formalism for performing joint analysis of cosmolog¬ 
ical datasets by introducing hyperparameter weighting 
terms, the yalues of which are determined directly from 
the statistical properties of the data. This approach is 
easily adapted to find the self-consistent uncertainties for 
data with ill-determined error -bars. Using Equations 29 
and 30 of lHobson et al. (2002) we calculate a likelihood 
function using the modified chi-squared statistic 


X 


2 


E 


(RMj - RM^od)" 


ln(2 7rcr(RM)2^t,i) . 


( 8 ) 

where RMi — RMmod is the difference between the i*** 
rotation measure and the model at that position. The 
second term inside the parentheses is required to cor¬ 
rectly normalise the likelihood and the total uncertainty 
on the RMi is giyen by 


cr(RM)?„t i = cr(RM)? -bexp[21n((5(RM)i)]. (9) 

Here we solye for ln((5(RM)) rather than directly for 
d(RM) so as to enforce positiyity in the scatter term 
(i.e., uncertainties cannot be negatiye). 


4. Results 

4.1. General comparison of model and data 

Before presenting the results of the MCMC analysis, 
it is useful to yisually compare the RM data shown in 
Fig. [3] with the simple shell models illustrated in Fig. 0 
Two key discriminators stand out in the behayiour of 
the models. Firstly, the difference in rotation measure 
between the interior and the peak of the shell (ARM, 
illustrated by the black line in Figures El is a strong 
function of the compression factor X, assuming other 
parameters are fixed. From Fig. [3]- right we see that the 
measured ARM in the profile of the northern Gum Neb¬ 
ula is ^ 350 rad m“^, restricting the compression factor 
to A < 4. Models with higher yalues of X result in much 
greater ARM for all reasonable yalues of Ue, /, Bq and 
0. Secondly, the longitudinal RM-gradient is directly 
related to the pitch angle of the magnetic field. This be¬ 
hayiour is due to the close proximity of the nebula, so 
that sight-lines from opposite sides are not parallel and 
intersect a uniform field at different angles. If the or¬ 
dered magnetic field is directed along the plane of the 
sky at the nebula’s centre (0 = 0°), we would expect to 
measure equal positiye and negatiye RMs on either side 
of the central longitude. For a magnetic field pointing di¬ 
rectly towards (0 = 90°) or away (0 = —90°) from the 
Sun the RMs would display symmetric positiye or nega¬ 
tiye patterns, respectiyely. We see mostly positiye RMs 
towards the Gum Nebula, with a slight positiye gradient 
towards lower Galactic longitudes (see Fig. m. From an 
examination of the grid of models shown in Figure El we 
can conseryatiyely state that the ordered magnetic field 


is pointing towards the Sun at an angle 0 > 20°. This 
is because the negatiye peak in RMs at positiye Galactic 
lon gitudes is absent fro m models with 0 > 20°, and from 
the Taylor et al. ( 20091) RMs. In ( 14.21 below we quantify 
these assertions using fits to the data. 

4.2. Fits to the model shell 

Here we present the results of fitting the mode l de¬ 
scribed in (13.2l to a subset of the Taylor et al. ( 20091 ) RM 
catalogue. The RM data included in the fit are outlined 
by the wedge-shaped box in Fig. E]- left and was selected 
to bracket the nebula aboye h > 5° (excluding the ‘stalk’ 
region and a handful of negatiye outliers - see (12.2.111 . We 
fixed the centre of the model to (Z, b) = (258.0°, —6.6°) 
so that the circumference of the shell corresponds to the 
sharp outer edge seen in the Ha data (see Fig. El- The 
mean of the marginalised likelihood distribution is a good 
estimator of the best fitting yalue for each parameter; 
these are reported in Tableland described below. 

We initially ran our MCMC fitting procedure assuming 
a flat, large-scale background of RMbg = —26.4radm“^ 
and with all other parameters free, except distance, 
which was fixed a.t D = 450 pc and electron density, 
for which a prior of rie = 1.4 ± 0.4 cm“^ was set (see 
( 13.311 . Likelihood distributions C and plots of RM are 
presented in Fig. El The triangular matrix of confidence 
contour plots illustrates how the free parameters inter¬ 
relate, while the histograms on the diagonal show the 
marginalised likelihood distributions for individual pa¬ 
rameters. Of particular note are the distributions for /, 
He, Bq, dr and X. The likelihood distribution for the fill¬ 
ing factor is very broad, only constraining / > 0.25, be¬ 
low which the marginalised distribution drops off rapidly. 
The curved and elongated confidence contours between 
Bq and / mean that these two parameters are highly 
degenerate, and that the strength of the magnetic field 
is not well determined in the absence of an indepen¬ 
dent estimate of /. The prior on electron density has 
the effect of constraining the likely range of rig values 
and eliminates most of the degeneracy between Ug and 
Bq. Gonfidence contours between dr and Pouter are el¬ 
liptical in shape, indicating that the RM-data do not 
pinpoint the radius and thickness of the shell indepen¬ 
dently. Nonetheless, the range of values for each parame¬ 
ter is small and their absolute values are well constrained. 
The marginalised C distribution for the compression fac¬ 
tor X exhibits a narrow profile centred on 1.1, imply¬ 
ing that the gas within the shell has not been signifi¬ 
cantly compressed. The hyperparameter characterising 
the additional scatter on the RMs is well determined at 
(5(RM) = 75.1i'j2’yradm“^, as shown by the Gaussian 
form of its marginalised likelihood distribution. By def¬ 
inition the hyperparameter procedure adjusts d(RM) so 
that = 1.0, thus <5(RM) can be though of as a proxy 
for x^ when comparing ‘goodness-of-fit’ between models. 

We ran the MGMG analysis again after correcting 
the data f o r the large- scale RM-gradient s mod elled by 
Sun et al. ( 20081) and Ijansson fc Farrar ( 2012h . The 
shape of the gradients is similar for both modes (illus- 





















































Table 2: Results of fitting the ionised shell model to the RM catalogue. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Parameter 

Symbol 

Unit 

Notes 


Assumed Background Level 





Flat 

Sun et al. 12008) 

■Tansson & Farrar 12012) 

Distance 

D 

pc 

Fixed 

450 

450 

450 

Background RM 

RMbg 

rad m~^ 

Fixed 

-26.4 

-64.0 

-36.6 

Filling factor 

/ 

- 

Free 

0 4 

-0.2 

-0.1 

-0.1 

Angular radius 

pouter 

deg. 

Free 

22 . 7+01 

22.7 to l 

22.7 tol 

Shell thickness 

dr 

pc 

Free 

20.2 tj® 

18.5 

18.51?:? 

Field angle 

0 

deg. 

Free 

55 til 

43 19 ® 

551?? 

Field strength 

Bo 


Free 

0 0 +6.1 

0-0 -4.0 

0 Q +4.9 
—2.2 

3.91?:? 

Compression factor 

X 

- 

Free 

1 1+0.5 
-0.3 

6 - 01 ?:? 

6 . 81 ?:? 

Electron density 

He 

cm“® 

Prior 

1.4 

1 0 +0.4 
-0.4 

1 0 +0.4 
-0.4 

Additional RM Scatter 

5(RM) 

rad m~^ 

Free 

75.112.? 

71.01?:? 

70.01?:? 


trated in Fig. a and removing these backgrounds has the 
effect of decreasing the RM-values towards lower Galac¬ 
tic latitudes. Due to the orientation of the selection box, 
this results in a decreased RM-signal towards the centre 
of the Gum Nebula compared to the edge. The parame¬ 
ters of the best-fitting models to the gradient-subtracted 
versions of the RM-catalogue are presented in columns 
( 6 ) and (7) of Tabled and illustrated in Figs. ITOl and fTTl 
The results of both MGMG fits are identical within the 


errors , so we refer to the version assuming the ISun et al 
( 2008l i background in the following discussion. 

Gomparing the results of fits to the flat- and gradient- 
subtracted data, the most significant difference is in the 
compression factor X, whose value changes from l.llg 3 
to respectively. This change is purely a result 

of the smaller difference between RMs in the interior 
of the nebula and the peak. X is much less well con¬ 
strained by the gradient-subtracted data as the difference 
between the interior and exterior levels approaches the 
scatter on the data. The confidence contours in Fig. [in] 
also show that X is more degenerate with / and Bq. 
The higher compression factor is balanced by correspond¬ 
ing small decreases in filling factor /, shell thickness dr, 
field strength Bq and electron density rig. The filling 
factor / is slightly better constrained, leading to a cor¬ 
respondingly more precise value for Bq = 

The uncertainties on fitted values of magnetic field angle 
0 are large (typically ^ 12 °), however, the fitted an¬ 
gles are broadly similar for all three fits. The value of 
(5(RM) is lower by 4.1 radm“^ in the gradient-subtracted 
fit, implying that the model is a better match to this 
data. However, the absolute change is equivalent to only 
~ 1.5cr between the two MGMG runs. We discuss the 
implications of the results in (J5l 


4-2.1. Consistency checks - pulsars 

Although fewer in number than extragalactic sources, 
pulsars with well determined distances are useful in 


checking the results of an extragalactic RM or EM anal¬ 
ysis. The interaction between free electrons and pho¬ 
tons introduces a differential time delay At across the 
observational bandwidth Az/. The delay is a function of 
frequency v and is characterised by the dispersion mea¬ 
sure according to {At/fis) = 8.3 Ai/(DM/cm“^pc) 
where all frequencies are in GHz. The measured DM of 
a pulsar is related to the electron density via 

npsr 

DM = / Uedl cm“^ pc. (10) 

J ohs 

Assuming the same volume filling factor / and path 
length L as before. Equation [TUI can be written as 

DM = 7 * 6 /A cm“^ pc. ( 11 ) 


With suitable observations and by combining Equa¬ 
tions m in and [TT] we can solve for /, B\\ or rig along the 
line-of-sight to a pulsar. For example, the average line- 
of-sight magnetic field strength is given by i?|| = g 
the electron density inside the clumps by rig = and 
the filling factor by / = 

The RMs, DMs and accurate distances (where avail¬ 
able) of selected pulsars towards the upper Gum re¬ 
gion have been presented in Fig. | 6 l A few general 
trends are worth noting: on average the DMs increase 
towards the Galac tic mid-plane as the electron density 
peaks at 6 « 0° ( Gaensler et al. 2008li . After taking 
the latitude dependence into account, the DMs of pul¬ 
sars towards the nebula appear to be enhanced when 
compared to sight-lines outside of its border. Typi¬ 
cally, pulsar sight-lines just inside shell edge have DMs 
of 100-150pccm“^ compared to 20-60pccm“^ out¬ 
side. The DM of the Vela pulsar (J0835—4510) is 
68 pccm“^, greater by 19pccm“^ compared to the pul¬ 
sar J0737—3039A, which lies just outside the nebula at 
distance of 1.1 kpc. The difference is representative of 
how much dispersion is created by one wall of the neb¬ 
ula. 
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Fi g. 9.— Main-figure : Triangular grid of confidence plots for the Gum Nebula, produced by fitting the ionised shell model to 
the lTavlor et al.l (120091 ') RM-catalogue assuming a flat background (also see Table[2]). The blue, red and green contours represent 
the Icr, 2a and Sn confidence intervals, respectively. The filled histograms on the diagonal show the marginalised likelihood 
distributions for each free parameter in the model. Inset: The rotati on measure image ( top) and azimuthal profile (bottom) for 
the best fitting model over-plotted by the selected RM data from the lTavlor et al.l (l2009l ) catalogue. Note the compressed range 
in angular offset on the x-axis of the inset panel, compared to the equivalent plots in Fig. [71 
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Fig. 10.— As for Fig. |9l but assuming a background as given bv ISun et al.l ll2008l ). 
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Fig. 11.— As for Fig. |9l but assuming a background as given bv ijansson fc Farrar! ll2012h . 


The clumpy electron density derived from rig = dm 
using only pulsars above b > 5° inside the nebula varies 
from 0.7cm“^ to 2.2 cm“^ with an average of 1.5 cm“^. 
Outside the nebula, but away from the plane, Uq (i.e., the 
ambient value) falls to values of 0.2 cm“^ to 0.9 cm“^. 
Given the large uncertainties these values are in agree¬ 
ment with our best-fitting models. 

Only a handful of pulsars inside the nebula have both 
RMs and DMs, allowing the determination of the average 
line-of-sight magnetic field strength i3||. Three pulsar 
sight-lines intersect the upper Gum region, away from 
the shell, and their DM values suggest they lie beyond 
the nebula. Values for derived from the pulsars range 
between 0.9/rG and 2.4/iG, significantly lower than the 
best-fit value to the flat-background data (S.SI^^oMG), 
but consistent with the value found when fitting the 
gradient-subtracted RM-data fiG). Some of the 


discrepancy may be explained by our choice of filling fac¬ 
tor, which is not well-constrained in any of the results. 
The fitted value of Bq is highly dependent on / and a 
values of / ~ 0.5 would bring our models into agree¬ 
ment with i?|| derived from pulsars. Determining the 
filling factor from a pulsar requires a good estimate of 
the path length and hence of the distance to the pulsar. 
Unfortunately, no pulsars with well measured distances 
lie towards northern part of the Gum Nebula, thus we 
do not attempt to estimate / at this time. 

In summary, we find that the pulsar data are consis¬ 
tent with our results, and especially favour the datasets 
which have had a model large-scale Galactic RM signa¬ 
ture subtracted (i.e., columns 6 and 7 in Table [21 fits in 
Figures HO] and . 

























































5. Discussion and further analysis 

The best-fitting shell models have some interesting im¬ 
plications, especially for the local direction of the ordered 
magnetic field and the fitted compression factor. We dis¬ 
cuss the results below, but start by noting the limitations 
of the model and the data. We further analyse the results 
by comparing the expected radio-continuum signature of 
the best-fitting models to the diffuse S-PASS 2.3 GHz 
data. 


5.1. Limitations of the simple shell model 

The simple ionised shell model presented here has a 
number of limitations and assumptions that should be 
considered when interpreting the results. We have al¬ 
ready noted in M.2I that the compression factor X is 
sensitive to differences in RM between the rim, exterior 
and interior of the nebula. The vertical orientation of 
the selection box and the relatively narrow portion of 
the nebula sampled by observations (a ~ 76° pie-shaped 
sector, see Fig. [3]- left) mean that latitudinal gradients 
in RM affect the fitted value of X in particular. Simi¬ 
larly, the fitted angle of the magnetic field 0 is directly 
dependent upon the observed longitudinal RM-gradient, 
as the ordered field runs parallel to the Galactic disk 


( Mathewson fc Fordlll970l : Han et aLll2006l) . Thus, iso¬ 
lating the RM-signal of the Gum Nebula is a critical step 
in our analysis and involves subtracting RMs due the 
large-scale Galactic background, and the smaller-scale 
magneto-ionic material along the line of sight. A residual 
RM-gradient remaining within the data would skew the 
values of X and 0 derived from our MGMG analysis. 

The lack of co verage below Dec. = —40° in the 
Taylor et al.l (200^ RM catalogue makes determining an 
accurate large-scale Galactic background difficult, so we 
have tried two approaches: subtracting a flat background 
and subtracting a model RM-gradient. Exterior to the 
Gum Nebula {b > 12°) the RM-data are small and neg¬ 
ative, consistent with a homogeneous background. How¬ 
ever, the RM values must increase towards the Galac¬ 
tic mid-plane, as the electron density is known to fall 
exponentially with increas ing latitude (jCordes &: Lazio 
2002 : Gaensler et al. 20081) . As discu ssed in ii2. 2. II pre ¬ 
viou sly' the Galactic RM-m odels of ISun et al.l (120081 ) 
and IJansson fc Farr an (120121) represent the best exist¬ 
ing estimate of the RM distribution due to the bulk 
of the Galaxy behind the nebula. Neither model is a 
good fit to the local RM distribution, poorly match¬ 
ing RM-structures on scales of ~ 10° in the vicinity of 
the Gum Nebula. H owev er, it is encouraging tha t both 
the Sun et al. ( 20081) and l.Tansson fc Farrail ( 2012h mod¬ 
els have similar gradients so we believe the large-scale 
morphology to be reliable, but not the local calibra¬ 
tion. Therefore, towards the mid-plane, RMs with the 
flat background subtracted constitute an upper limit on 
signal from the Gum Nebula. 

On small scales the division of RMs into ‘backgronnd’, 
‘Gum Nebula’ and ‘other object’ categories is necessary 
to obtain a clear RM-signal (see (12. 2. H and Fig. |3]). This 


identification procedure draws on all of the available data 
to make informed decisions, but the process is still some¬ 
what subjective. Residual RMs from unidentified dis¬ 
crete objects may still be present in the data, or the 
identified objects may extend behind the footprint of the 
Gum Nebula. For example, the RM signature of a small 
Hii region (< 2°) overlaid on the rim could be erro¬ 
neously fitted as a gradient in I or &, leading to a sys¬ 
tematic errors in 0 or A, respectively. Future surveys 
that deliver a more accurate and densely sampled grid 
of RMs covering the southern sky are required to resolve 
remaining ambiguities. 

It is clear from the filaments visible in the Ha map 
that Ue is structured on scales down to the 6' resolution 
of the image. This clumpy distribution of electrons is 
acconnted for in the model using a global volume filling 
factor /, leading to an occupation length / L for all sight¬ 
lines. Variations in rig on scales much smaller than the 
beam lead to fluctuations in RM, which manifest as an 
additional uncertainty on the RMs, codified as d(RM) 
in the model. Because the spatial sampling is coarse 
(~ 1/degree^), the value for d(RM) = 75.lt2'7radm”^ 
is an upper limit on the true scatter in rotation measure. 
The high value may also reflect genuine scatter of ISM 
properties between the model and the data. 

We have assumed the electron density profile ne(r) 
within the shell is constant as a function of radius. This 
is in line with the d escription of a wind-blown bubble (see 
IWeaver et al. 1977 . Fig. 3), or with the physics o f an ex¬ 
pand ing ionisation front in an evolved H ii region ( Drainel 
l201l[) . However, a constant ne(r) is inconsistent with 
the density profile of a Sedov-phase SNR , which increases 
from the centre towards the shock-front (|van der Swaluw 
200 ih . The measured density profile of the Gum Nebula 


presented in Fig. [8] agrees with a constant ne(r) within 
the errors. The density is slightly enhanced towards the 
front edge of the shell, but only at a ~ Icr level, so we 
consider this assumption reasonable. 

The model accounts for compression at the edge of the 
shell using a factor X by which both the ambient density 
no and the component of the magnetic field tangent to 
the shell surface {B±) are amplified. The model assumes 
that the gas inside the shell is 100 percent ionised by the 
powering source or the passing shock. The magnetic field 
component that produces the RM signature (73||) is given 
by the projection of B± and the radial component. Bn, 
onto the line-of-sight. As a computational convenience 
i?l| is assumed to have a constant strength throughout 
the thickness of each shell wall; an assumption which is 
valid only for a thin shell. For the Gum Nebula, the 
ratio R/dr « 15 and the best-fit compression factor is 
low {X < 10), so this assumption is acceptable. 

The model has implicitly assumed that the lines of the 
ordered Galactic magnetic field are parallel to each other 
and to the disk of the Galaxy. If they loop, converge or 
diverge significantly within the nebula (260 pc) then a 
much more sophisticated treatment is required, coupled 
with a more finely-sampled grid of rotation measures. 
Such an analysis is beyond the scope of this paper but 



























































should be considered with future datasets. 


5.2. Orientation and strength of the ordered 
Galactic magnetic field 

When viewed face-on to the Galactic disk, the direc¬ 
tion of the large-scale Galactic magnetic field is charac¬ 
terised by the pitch angle, defined as the deviation from 
a circular path around the Galactic centre and given by 
p = tan“^(i?rad/i3az), where Bradiai and Baz are the ra¬ 
dial and azimuthal components of the ordered field, re¬ 
spectively. In external spiral galaxies the magnetic field 
lines are observed to closely follow the spiral arm pattern, 
but the field strength is o ften greatest in the inter-arm re¬ 
gion (see the examples of B^eck et al.ll2005l Fletcher et al 


2 OO 4 I and IPatrikeev et all 200fili . In the Milky Way, the 


ordered disk field is directed parallel to the disk with a 


X (kpc) 

-10 -50 5 10 



Fig. 12.— Top panel: Cartoon showing the location of the 
Gum Nebula in the disk of the Milky Way. The coloured 
lines trace the pol ynomial-logarithmic spiral arm models of 
IHou fc HanI l|2014h . A dashed magenta line plots the solar 
circle at a Galactic radius of 8.5 kpc. Bottom panel: Illustra¬ 
tion of the fitted magnetic field orientation around the Gum 
Nebula. The plot covers the region inside the grey box in the 
upper panel and the Gum nebula is shown by a black circle 
(to scale). The three thick arrows piercing the Gum Nebula 
show the median (black) and ± 1-sigma ranges (grey) of field 
orientation found by fitting the RMs in 1)4.21 The field of thin 
grey vectors depicts the orientation of a spiral magnetic field 
with a pitch-angle of p = —7.2°. 


typical strength of Bq = 1.5 — 2^G ( Han et al. 20061) . 
The pitch angle of the field in the disk has been estimated 
by multiple authors using a variety of techniques and has 
been found to lie between —6° and —11.5°, depending on 
the method used and the volume of the Galaxy observed. 
There is also some evidence that p may have a radial 
dependence, decreasing to almost zero at galactocentric 


radii greater than the solar orbit (IVan Eck et al.M2011 


Jansson fc Farraill2012h . Table [3] summarises the results 


of individual studies in the literature. 

The pitch angle p. Galactic longitude I and fitted mag¬ 
netic field angle 0 are related by simple geometry via 
I — 180° = 0 -I- p. At the Galactic longitude of the Gum 
nebula {I ~ 258.0°) the median pitch angle of p ~ —7.2° 
from Table [ 3 ] implies an ordered field pointing almost di¬ 
rectly towards the observer (0 = 85.2°). Our best-fitting 
shell models presented in ij4.2l return field directions be¬ 
tween -1-43° < 0 < -1-55°, equivalent to a pitch angle 
range -1-23° < p < 4-35°, substantially different to pre¬ 
vious results from the literature. Taking the Tier limits 
for all models, the local pitch angle is constrained by 
our data to 4-7° < p < 4-44°. This range represents 
the pitch angle of the uniform ambient field local to the 
Gum Nebula. Our results are illustrated in Fig. [121 which 
shows the position of th e Gum Nebula wi th respect to 
the Galactic spiral arms ( Hou fc HanI 20l3 ) and an ideal 
uniform field with p = —7.2°. All of the investigations 
referenced in Table [3] calculated the ‘global’ pitch angle 
averaged along the line-of-sig ht to pulsars, radio -galaxies 
or stars. All but one study ( Pavel et al.ll2012l ) covered 
a broad swathe in Galactic longitude. As such, the de¬ 
rived pitch angle is an average over a large fraction of the 
Galactic disk. By contrast, the method presented here 
probes only the magnetic field around the Gum Nebula, 
on scales of ~ 260 pc (the diameter of the nebula). 

Only a handful of other studies have used bub- 
bles as probes o f Ga lactic magnetic field structure: 
iKothes fc Browii ( 20091) studied two SNR in the Cana¬ 
dian Galactic Plane Survey ([Tavlor et al.l [2003 ) and 


found support for an azimuthal disk field, while [Ransom et al 


( 201(^ 1 performed a similar study using old planetary 
nebulae, but d erived only the li ne-of- sight field strength 
Most recently Whiting et al.l (20^ and ISavage et al.1 
( 2013h derived the angle of the magnetic field by mod¬ 
elling Hii region shells and found field directions com¬ 
patible with the mean Galactic field. 

We caution that some of the deviation in p may be due 
to systematic errors in the data. As discussed in 05.11 0 
is sensitive to longitudinal gradients in RM introduced by 
contaminating magneto-ionic material, or by the Galaxy 
in the background. We have taken all reasonable steps 
to identify and eliminate such contamination, however, a 
definitive correction requires much more finely sampled 
and accurate grid of RMs. 

The best fitting ambient magnetic field strength of 
Bq = 3.9l2'2 mG is within the range of 2 uG to 15 uG ob¬ 


served by similar studies of H i i regions ([Gaensler et al 


2001 : Harvev-Smith et al.ll2011 ). As shown in Fig.Hni B( 


is correlated with /, which is very poorly constrained by 


















































































Table 3: Studies of ordered magnetic field pitch angle in the literature. 


Pitch Angle 

Reference 

Notes 

-hl6° ±4° 

Inoue & Tabara 11981) 

Radio-galaxies, < 2 kpc, Orion arm 

-6° 

Vallee 11988) 

Pulsars, few kpc, Sagittarius and Perseus arms 

-8.2° ±0.5° 

Han & Oiao 11994) 

Pulsars (thin disk) and radio-galaxies (thick disk) ^ 3 kpc 

-8° 

Han et al. 11999) 

Pulsar RMs, ~ 15 kpc 

-7.2° ±4.1° 

Heiles 11996) 

Starlight polarisation, few kpc 

-11.5° 

Van Eck et al. 12011) 

Radio-galaxies, Galactic sector average 

-6° ± 2° 

Pavel et al. 12012) 

Radio-galaxies, average along 1 = 150° 


the data. The value of So above is reported for / = 0.3, 
the mean of the marginalised likelihood distribution. If 
instead / is set to the most likely value of / = 0.24 then 
Bq Ki b iiG. The strength of the field within the shell 
depends on the position, and varies between the ambient 
level and a maximum value oi X x Sq ~ 23 fiG when the 
Bq lies parallel to the edge of the nebula. 

In summary, the pitch angle of the ordered magnetic 
field threading the Gum Nebula (7° p ^ 44°) is 
significantly different to previous measurements, most 
of which were averaged over kiloparsec-sized volumes. 
Few small scale measurements of the field in the dif¬ 
fuse ionised medium exist, so this result may represent 
typical deviations on scale s of a few hundred parsecs. 
Indeed, iFrisch et al.l (12012!) measured even larger devia¬ 
tions in the ordered magnetic field in the vicinity of the 
Sun (< 40pc), consistent with a scenario where the local 
ISM is a fragment of the Loop I superbubble. Such de¬ 
viations have al s o bee n observed in external galaxies, for 
example HealdI ( 20121) detected a significant RM gradi¬ 
ent in the spiral galaxy NGC 6946, tracing a irregularity 
in the vertical component of the ordered magnetic field. 
The deviation is directly associated with a hole in the 
HI image and may be ubiquitous feature of star-forming 
galaxies. Expanding bubbles in the disk may also be re¬ 
sponsible for carrying the small-scale turbulent magnetic 
field into the halo, preventing quenching of the dynamo 
process and allowing the mean magnetic field to satu¬ 
rate at a strength compar able to equipartition w ith the 
turbulent kinetic energy ( Shukurov et al. 20061 ). More 
accurate and better-sampled RMs are required to con¬ 
firm our result and eliminate systematic uncertainties. 
The strength of the ambient field around the Gum Neb¬ 
ula is comparable to averag e values of 2 — 4 u G measured 
for the Galaxy as a whole (jHan et al.l 120061) and within 
H II regions. 


5.3. Implications of the fitted compression fac¬ 
tor 

The best fitting models presented in M.2l constrain the 
jump in density at the edge of the nebula, assuming the 
shell is 100 percent ionised (by stellar radiation in the 
case of a H ii region or wind-blown-bubble, or by the 
shock-front in the case of a SNR). The fitted value for the 
compression factor assuming a flat Galactic background 
is X = 1.1 to ;3 and assuming a gradient is X = 6.0^ 2 ^. 
At the very least, both values imply that the gas within 


the shell is only moderately compressed compared to the 
ISM external to the nebula. 

The current consensus in the literature is that the 
nebula is an old supernova remnant (see )!l.ll). SNR 
pass through three distinct evolutionary phases ( Woltieij 


1972h before dissipating: 1) free expansion {t < 300 yr), 


where the swept-up mass is much less than the ejected 
mass and the expansion is dominated by the explosion; 
2) the Taylor-Sedov phase (300yr < t < 20000yr), where 
the swept-up mass dominates and the blast wave expands 
adiabatically and 3) the snow-plough phase (20000 yr < 
t < IMyr), when thermal cooling has become effective 
and the shock front decelerates, sweeping u p a dense 
shell. According to the widely used models of Chevalieil 
( 1974h the 260 pc diameter of the Gum Nebula implies 
an age of 1 Myr, which would be old indeed for an 
SNR. At this late stage of evolution the shock front is ex- 



SNR expansion is expected to slow down to the ambient 
sound speed (typically ~ 10kms“^) and merge with the 
ISM, although the exact details of this process are not 
clear ( Pittard et al.l 20031) . The expansion velocity mea¬ 
sured from optical spectroscopy towards the Gum Nebula 
is V S. 10kms“^ ( Srinivasan et al ] I1987L ISahu fc Sahul 
1993h . This slow speed and moderately low compression 


factor, combined with the flat rig profile derived in 
make it unlikely that the nebula seen in H a stems solely 
from a supernova origin. Instead, the detected H a emis¬ 
sion and RM-signature are consistent with an ionisation 
front moving at subsonic speeds into the ISM. Both a 
classical H ii region and wind-blown-bubble are bounded 
by an ionisation front, so we consider these models in 
turn below. We cannot completely rule out the old SNR 
origin as during the dissipation stage the shell likely re¬ 
expands, leading to a decrease in density and hence com¬ 
pression factor. However, at this stage, we would also 
expect the shell to loose cohesion as it merges with the 
ISM and this is not seen in the Ha data towards the 
Gum Nebula. 

The fundamental theory of e xpanding, o v er-pre ss urised 
Hii r egions was s et ou t by IStromgr^ ( 1939f ). Kahn 
( 1954h and Oorfj ( 1954h . After the initial formation 
phase the H ii region expands approximately isothermally 








































































( Dvson et alJ[l9^ 1 until pressure equilibrium is reached. 
For a high-mass 0-type star this does not happen within 
its stellar lifetime. The D-type ionisation front moves 
into the ISM at the local sound speed < 10kms“^ and 
has a three part structure, consisting of a thin layer of 
shocked neutral gas separating the unshocked neutral gas 
from the ionised gas within the H ii region. This classical 
description produces a spherical ionised region of approx¬ 
imately constant electron-density, at odds with the ob¬ 
served shell-like structure of the Gum Nebula seen in H a 
emission. One possible way to reconcile the model with 
the data is if the ionisation front is expanding into the 
walls of a pr e-existing cavity. T his explanation was first 
proposed by Revnoldj (19763) whereby the ultraviolet 
flux from the central stars is ionising the walls of a void 
and illuminating the ambient magnetic field in the ISM 
local to the Gum region. Such a cavity could have been 
formed by an older supernova explosion, or evacuated by 
an older generation of stars. 

The ionised cavity explanation is somewhat contrived, 
especially since the wind-blown-bubble model can nat¬ 
urally account for the structure and physical parame¬ 
ters of the Gum Nebula measured from observations to 
date. The star Puppis is known to drive a power¬ 
ful stellar wind, as is 7 ^ Velorum and the Vela OB2- 
associatio n (which may also l ie inside the Gum Neb¬ 
ula). The IWeaver et al. ( 1977 1 description of a bubble 
blown by a high-mass star predicts a shell-like H ii region 
surrounding a regio n of shocked stellar wind. Indeed, 
Weaver et al ] (ll977ll offer the Gum Nebula as a prototype 
wind-blown-bubble powered by the strong stellar winds 
from C Puppis. Figure 3 in their paper illustrates the 
temperature and density profile of a typical bubble. The 
density jump across the outer boundary of the H ii shell 
is X ss 1.5, broadly consistent with our results. In addi¬ 
tion, the density in the interior is low at Ue 0.05 cm“^ 
(as is known to be the case for the Gum Nebula) and in 
the shell is rie ~ 2.5 cm“^, com parable to me a surem ents 
in this work. Equation 69 in IWeaver et al.l (Il977ll de¬ 
scribes the density in the shell compared to the ambient 
upstream density 


= no{Vi + C^o)/Cl 


( 12 ) 


where V2 is the shell expansion velocity, Cq is the speed of 
sound in the ISM and Cg is t he sound speed in th e shell. 
For typical values used in the I Weaver et al. (19^ model 
of the Gum Nebula (V 2 ~ 8 kms“^, Cq ~ lkms“^, Cg ~ 
lOkms”^) then X = rie/no = 6.5 in line with our best 
fitting shell model assuming a background RM-gradient 
(see 112 . 2 . 11 ) . 

In summary, we believe our results point to the wind¬ 
blown-bubble model as the most likely explanation for 
the origin of the Gum Nebula. 


5.4. Pressure and evolutionary state 

The ratio /3th = Tth/Tmag of thermal to magnetic pres¬ 
sures in an Hu region can indicate whether the object 
is still evolving or has reached an equilibrium state. In 
a young H ii region, the thermal pressure dominates and 


drives the expansion of the warm ionised gas into the 
ISM, sweeping up ambient gas before it. If the region is 
threaded by a uniform magnetic field, flux-freezing in the 
ionised gas will resist expansion perpendicular to the field 
lines. Over time, as the Hii region expands, the thermal 
pressure decreases and the magnetic pressure increases, 
so ratios closer to unity imply an older object. Magnetic 
pressure is given by 

Pmag = .Bg/(S tt) dynecm”^, (13) 

where Bq is the total magnetic field strength in /rG. The 
thermal pressure in an ionised gas at temperature Te 
(in K) and density Ue (in cm“^) is 

Pth = 2nefcTe dynecm”^, (14) 


where k is Boltzmann’s constant. Assuming Te = 
8000 K and rie = 1.4 cm“^ for the Gum Nebula, we 
find Pth = 2.9 x 10“^^ dynecm”^ compared to Pmag = 
6.1 X 10“^^ dynecm“^. The ratio of thermal to magnetic 
pressure /3th = Pth/Pmag = 4.8 suggests that the dy¬ 
namics of the shell are dominated by thermal motions, 
i.e., the magnetic field is too weak to shape the overall 
morphology of the ionised gas. This result is in keep- 
in g with the magneto-hyd rody namic (MHD) simul ations 
of iKrumholz et ah (2003) and I Arthur et ah ( 201 lll . who 
find that the thermal pressure of the ionised gas shapes 
the evolution on time-scales of several Myr. Both simu¬ 
lations calculate similar field-strengths and gas-densities 
to what we derive for the Gum Nebula. Our values for 
Uf, , Bn and also sit i n the middle of the range found 
by iHarvev-Smith et al. ( 2011 ) in their survey of high- 
latitude evolved Hii regions. 

5.5. Radio Spectral Index 

The Gum Nebula is a prominent foreground feature 
in the recently released Planck all-sky ra dio-continuum 
maps ( Planck Gollaboration et al. l2014all . The 36° 
angular diameter shell is visible in the 28.4 GHz image 
at the lOcr level and in the 44.1GHz image at the 

3(7 level, implying a flat spectral index and a significant 
thermal component to the emission. A more quantative 
method of determining the spectral index of diffuse emis¬ 
sion is pr ovided by the tempe rature-temperature (TT) 
plot fe.g.. iTian & Leah3l20f)6ll . An x-y plot of the flux 
densities within a sampling aperture results in a scatter 
plot and the slope of a straight line fit to the data gives 
the spectral index S = Sov°‘. The main advantage of 
a TT-plot is that large scale offsets in the background 
emission are automatically compensated for, assuming 
that the background does not vary significantly within 
the sampling box. We find that the spectral index of the 
Gum Nebula shell at the brightest region is a = 0.2±0.2, 
based on a comparison of the 28.4 GHz Planck data and 
2.3 GHz S-PASS data; consistent with thermal free-free 
radio emission. We note that parts of the Gum Nebula 
are als o faintly visible in th e 408 MHz radio-continuum 
map of Haslam et al.l ( I982h , implying a mixture of ther¬ 
mal and synchrotron emission in places. 
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Fig. 13.— Top panel: Map of the 2.3 GHz polarised intensity for the upper Gum Nebula. Green rectangles show where 
profiles have been extracted, in Fig. and dotted lines outline the Gum Nebula and Antlia SNR. The profiles run from 
(/,6) = (268.0°, 17.7°) to , (263.1°, 11.9°) and (259.1°, 18.9°) to (259.1°, 11.7°) for profiles 1 and 2, respectively. Bottom panel: 
Map of the dispersion in polarisation angle The value of cr{ip) in each pixel has been calculated from the standard deviation 

within a beam-sized aperture. Dark canals correspond to regions where the polarisation angle varies by close to 90° within a 
beam. The red bars over-plotted on the image every 7th pixel illustrate the orientation of the polarisation vectors (magnetic 
vector, not corrected for RM). 
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Fig. 14.— Two profiles extracted from the map of the 
2.3 GHz polarised intensity. The tracks along which the pro¬ 
files have been extracted are illustrated by the green boxes in 
Fig. [island the starting positions marked by circles. 



Fig. 15.— The depolarising effect of the spatial gradient in 
RM at the inner and outer edge of an ionised shell. The dot¬ 
ted/red line shows the gradient in rotation measure (dRM/dr) 
in units of rad m~^ beam” as a function of position, cross¬ 
ing the shell from outside to the interior. The solid/green 
line shows the corresponding degree of polarisation calculated 
from Equation 1171 


Further investigation of the spectral index through¬ 
out the Gum Nebula is beyond the scope of this work 
and should be the subject of a separate paper. There 
now exist many wide-angle maps of radio-continuum 
emission covering the Gu m Nebula, includi ng histor- 
ical data (e.g., 4 5 M Hz: Maeda et al.l 19^ 1.4 GHz: 
Reich et al.l 2001 and Galabretta et al. 20l4 300 MHz - 


1.8GHz: Wolleben et al. 2009 . 5 GHz: King et al. 20101 
23-94 GHz: Bennett et al. 2013ll and new diffuse polar¬ 
isation maps from Planck and the Murchison Widefield 
Array ( Tingav et al. 2013ll . Future investigations com¬ 
bining these datasets will be capable of disentangling 
emission due to synchrotron, free-free and ‘spinning dust’ 
processes across the region. 

5.6. Polarised 2.3 GHz radio-continuum emis¬ 
sion 

The properties of the diffuse polarised 2.3 GHz radio 
emission provide complementary information to the RMs 
of background radio-galaxies. In particular, analysis of 
the polarised intensity and the angle of the linear polari¬ 
sation vector can yield information on the geometry and 
the level of turbulence in the ionised gas. 

Fig. Ill] presents two views of the polarised 2.3 GHz 
emission centred on the upper shell of the Gum Neb¬ 
ula. The top panel displays a high resolution image of 
the polarised intensity (P = + The rim of 

the Gum Nebula stands out as a broad 2° wide) 
band of depolarisation across the centre of the image. 
In addition to the Gum Nebula, two other objects have 
been identifie d in the field. The ed ge of the Antlia 
SNR ^2.2.11 ( McCullough et al. 2002ll is visible in the 
upper-left quadrant as an arc of weaker depolarisation 
and narrow canal-like features. Such canals trace re¬ 
gions where the polarisation angle ■i/' varies significantly 
across a telescope beam, leading to d epolarisation in the 
receiver ( Fletcher fc Shukurov 200^. The Antlia SNR 
overlaps the Gum Nebula between 163“’ < I < 172“’, 
where the Ha emission is brightest (see Fig. [Ij. Also 
visible are prominent depolarisation canals from a pair 
of shells in the upper-right qu adrant of the image (see 
il2.2.1l) . lacobelli et al.l ( 2014ll have analysed t he spa¬ 
tial polarisation gradient (jGaensler et al.l 1201 ill in the 
S-PASS data and have identi fied these features as the 
signature of weak shocks (see iBurkhart et al.l 12012 ). It 
is not known if either object is physically interacting 
with the Gum Nebula, or is simply seen in projection 
along the line-of-sight. The bottom panel plots the po¬ 
larisation angle of the electric vector (red bars) over an 
image of the spatial dispersion in the polari sation an¬ 
gle (cr(^) = stdevl^Vjl. iHildebrand et aklboOfll l. Maps of 
aitl;) have been shown to highlight depolarisation canals 
( Planck Collaboration et al.l 120145) and are useful way 
of visualising where the polarisation vectors are homoge¬ 
neous or heterogeneous on the sky. Within the Gum Neb¬ 
ula and Antlia SNR the polarised intensity is patchy and 
the polarisation angles are chaotic in comparison to the 
slowly varying distribution of angles outside their bor¬ 
ders. We interpret this as a ‘scrambling’ of the smooth 



















































































synchrotron emission from the Galaxy in the background 
by Faraday screens associated with each object (e.g., 
ICarretti et al. 2ni3b). It is clear from preyious studies 
Duncan et al Ifl^ and radio-continuum data (; 


{e.g.. liJuncan et al ana raaio-continuum data (see 

s that emission from the Gum Nebula is dominated 
by thermal processes. The shell of the nebula does not 
emit significant amounts of synchrotron radiation and so 
can be analysed as a pure Faraday screen. Assuming a 
smooth synchrotron background of polarised synchrotron 
radiation from the Galaxy, we can quantify the effects of 
the screen by examining how the leyel of polarisation 
changes across the edge of the Gum Nebula. 

Fig. [14] presents two polarised profiles extracted from 
the polarised intensity map. Prohle 1 cuts across the 
brightest region of shell, where it oyerlaps the weaker 
Antlia SNR and Profile 2 has been extracted from 
the upper part of the shell. Profile 1 drops from a 
high of 55mJybeam“^ outside the shell to a low 
of 15mJybeam“^ inside the depolarised region. The 
typical root-mean-squared intensity in the P image is 
2.2mJybeam“^ and the broad depolarisation band ex¬ 
hibits polarised emission at a ~ 5cr leyel. Assuming the 
outer yalue represents the intrinsic polarised intensity 
Pi, then the degree of depolarisation implied by the drop 
to Pq = 15mJybeam“^ is p = Po/Pi ~ 0.27. The drop 
in intensity and degree of depolarisation is similar for 
Profile 2. 


5.6.1. Depolarisation 

The c auses of dep ol arisation haye bee n described in 
detail by [Burn ( 1966h . Tribble! (1991) and Sokoloff et al 


(11991^ . The root cause in all cases is cancellation be¬ 
tween polarisation yectors oyer some ayeraging interyal 
in time, space or frequency. 

Bandwidth depolarisation occurs when Faraday rota¬ 
tion causes the polarisation angle to yary across a fre¬ 
quency ayeraging window Av. The degree of depolarisa¬ 
tion due to frequency ayeraging is 


P = 


sin Atjj 


Alp 


(15) 


where the change in angle across a band centred on vq 
is giyen by Aip = —2RMc^ /z^g. In the Gum Neb¬ 
ula the maximum RM detected is ^ 350 radm“^ so the 
expected angle change oyer the 244 MHz bandwidth is 
Alp = 72° and the resultant depolarisation is negligible 
at p = 0.75. A RM of 870 radrn”^ would be necessary to 
completely depolarise S-PASS data. 

The most likely depolarisation mechanism affecting the 
S-PASS data is beam depolarisation. This is caused by 
yariations in P|| or rig on scales much smaller than the 
beam, s cattering the polarisation angles on adjacent lines 
of sight. [Burn ( 1966!) quantified this effect in the simplest 
case of a uniform slab and found 


p = exp(-2cr^MA‘‘)) (16) 

where ctrm is the RM scatter within a beam after mea¬ 
surement errors haye been accounted for. If small-scale 


random fluctuations are solely responsible for the ob- 
seryed depolarisation (p = 0.27) then Equation |T6| pre¬ 
dicts an excess scatter of ctrm = 47radm“^. In il4.2l 
we found that the best-fitting model implied an addi¬ 
tional scatter of ctrm = 78.6radm“^ (called (5(RM) in 
Table [5]). This fitted yalue is an upper-limit on ctrm 
as the RM sampling grid is yery coarse at ^ 1/degree^, 
compared to the beam FWHM of ©beam = 10.75'. We 
can conclude that the data is at least consistent with a 
large fraction of the depolarisation being due to random 
fluctuations in i3|| or Ug. 

While the ayerage drop in P can be explained (at 
least in part) by fluctuations within an ionised Faraday 
screen, the shell of the nebula also contains depolari¬ 
sation canals. These are typically one beam in width, 
close to 100 percent depolarised and tend to be aligned 
pa rallel to the edg e of th e nebula. First discoyered 
by Hayerkorn et al. ( 2000!) . seyeral authors in the last 
decade haye studied origin of suc h canals and explored 
their use as a diagnostic to ol (' e.g.. iFlet cher fc Shukuroy 
20061 Gaensler et al.l 2011 a.nd Burkhart et al.ll2012l l. In 
particular, iGaensler et al.l (|201lh calculated the spatial 


gradient of the complex Stokes yector = (^, l}) 
whose magnitude |VP| describes the rate at which the 
polarisation yector traces out a path in the Q-U plane 
when moying along a spatial track at a constant rate. 
|V P| is inyariant under arbitrary rotations or transla¬ 
tions (unlike P or 'tp) and images of |V P| reyeal a net¬ 
work of filaments in the ionised gas (see llacobelli et al 


I 2 OI 4 I for the |VP| of the S-PASS data). In a pure 
Faraday screen these filaments haye been shown to 
trace spatial cusps or jumps in rig or most likely 
caused by shock-front s or turbulent motions in the gas 
([Burkhart et al.ll2012[) . Depolarisation canals like those 
in Figure [131 are a subset of filaments that cross the origin 
in the Q-U plane. The greatest concentration of canals 
occur within the rim of the nebula, lending weight to our 
conclusion that turbulent fluctuations in Ue or i3|| are 
responsible for the depolarisation. 

The most prominent canals run along the inner and 
outer edges of the depolarised rim and can be explained 
by the intrinsic RM-gradient at the edges of the ionised 
shell. The amount o f depolarisation prod uced by an RM- 
gradient is giyen by ISokoloff et al.l ( 19981) as 


Pg = exp 


1 /dRMV 
ln2 \ dr ) 




(17) 


assuming a Gaussian beam which resolyes the gradient. 
Fig.[T5|plots the RM-gradient and the depolarisation fac¬ 
tor calculated from the Equation flTl for a profile cross¬ 
ing towards the interior of the ionised shell. From the 
plot we see that depolarisation only becomes significant 
(p < 0.6) close to peaks in the gradient. The equation 
breaks down for resolyed gradients, howeyer, it is clear 
that narrow depolarisation canals are predicted at the 
leading and inner edges of the shell. 

In conclusion, the polarisation and depolarisation 
properties of the 2.3 GHz S-PASS data are in keeping 












































































with the simple ionised shell model put forward in 
and support our assertion that the Gum Nebula is acting 
as a Faraday screen. 

5.7. Comparison to previous studies 

The first dedicated magnetic field measurements of 
the G um Nebula were obtained by IVallee fc Bignell 
(1983) via linear polarisation observations of 35 back¬ 
ground extragalactic radio sources. Prior to that work, 
large-scale rotation measure excesses in the area were 
attributed to a tangential view of the local Orion - 
spur spiral arm ( Simard-Normandin &: Kronberg 1980l) . 


Vallee fc Bismelll (1983) claimed that the distribution of 
RMs on the sky were not consistent with the arm model 
but were a good ma t ch for the old SNR model first pre¬ 
sented by Revnoldj (1976b). Their derived line-of-sight 
magnetic field strength of ^ 1.3/rG suggested that a 
‘snow-plough’ effect alone was responsible for sweeping 
up gas, and hence th e magnetic field lines. In later work, 
Duncan et al.l (1993) cast doubt on the significance of 
the Vallee model, pointing out that the statistical uncer¬ 
tainty in the data used therein was comparable to the 
mean RM value. The mode l presented in th i s pap er is 
broadly consistent with the Vallee fc Bignell (1983) re¬ 
sult, but is considerably more sophisticated and includes 
much better sampled measurements of RM and Ue- We 
also derive independent values for the shell th i ckness 
and compression factor, which IVallee &: Bignelll (|l983t l 
did not provide. 

Magnetic field strengths in ionised bubbles have been 
meas ured by a number o f rec ent studies in the l itera- 
ture. Whiting et al.l ( 20091) and Savage et al.l ( 20131) used 
a similar technique to the one presented here to study the 
bubble surrounding the Cygn us OBI associa t ion an d the 
Rosette nebula, respectively. IWhiting et al.l ( 20091) sug¬ 
gested that the observed Faraday ‘anomaly’ was caused 
by a wind-blown bubble, but with only nine RMs they 
could not confirm the compression factor predicted by 
the st rong shock. On the other hand, Sayage et al.l 
( 2013[) modelled RMs seen through the Rosette nebula 
as a limb-brightened ionised shell and obtained a con¬ 
siderably better fit when fixing X = 4 compared to 
X = 1. Recently, iHaryey- Smith et HI ( 2011 ) studied 
the line of sight magnetic fields in fiye large-diameter 
H II regions offset from the Galactic plane. They deriyed 
field strengths from ^3 — 11 /rG, but found no eyidence 
of compression at the edges of these relatiyely eyolyed 
H II regions. 

6. Summary and Conclusions 

We haye deyeloped a simple model of the Gum Neb¬ 
ula as an expanding ionised shell threaded by a uni- 
fo rm magnetic f ield. Drawing upon the R M catalogue 
of Taylor et al.l ( 2009t) and the H a image of iFinkbeined 
(l2n03^ . we used a maximum-likelihood MCMC analysis 
to deriye the magneto-ionic shell parameters in the upper 
hemisphere of the nebula. We compared the best-fitting 
models to polarised 2.3 GHz radio-continuum emission 


from the S-PASS project. Our conclusions are as follows: 

1. The RM and EM data coyering the upper hemi¬ 

sphere of the Gum Nebula (6 > 5‘’) are well-fitted 
by a simple ionised shell. Assuming a lar ge-scale 
RM-background from the lSun et al.l ( 20o3 ) model 
of the Galaxy, the best-fitting shell has an angu¬ 
lar radius (/<outer = shell thickness dr = 

18.5ti;4PC, ambient magnetic field strength Bq = 

electron density rie = l.d+o^cm”^ and 
filling factor / = 0.3lo;^ 

2. We constrain the pitch angle of the uniform mag¬ 
netic field to yalues oyer the range -1-7'’ < p < 
-1-44'’, significantly different from preyiously de¬ 
riyed yalues (p ~ —7°) ayeraged oyer much larger 
yolumes of the Galactic disk (scales of seyeral kpc 
yersus 260pc for the Gum Nebula). Our fitted 
yalues are sensitiye to contamination of the RMs 
by interyening magneto-ionic objects, howeyer, we 
haye corrected the catalogues to the full extent al¬ 
lowed by the ayailable data. This represents one of 
the few measurements of local magnetic field ori¬ 
entation in the Milky Way. 

3. We find that the compression factor X = Ue/rio 
at the edge of the Ha shell is X = 6.0as - 
suming an RM-background from Sun et al. ( 20081) . 
This yalue is much lower than expected if the Gum 
Nebula were an old SNR cooling radiatiyely. We 
belieye that the most likely explanation for the 
nebula is a wind-blown-bubble driyen by a cluster 
of high-mass stars. The slow expansion yelocity 
(< 10kms“^), low excitation conditions and lack 
of radio-synchrotron emission from the rim is con¬ 
sistent with our hypothesis. 

4. The strength of the ordered magnetic field Bq is 
not well measured as it is degenerate with the ill- 
constrained filling factor /. We deriye a value of 
Bq = 3.9j'j2’2 mG, in line with the strength of the 
ambient Galactic field and also comparable with 
values measured towards Hii regions by previous 
authors. 

5. Viewed in 2.3 GHz radio-continuum, the upper 
shell of the Gum Nebula exhibits a distinctive band 
of depolarised emission. We find that the dominant 
depolarising mechanism is likely due to fluctuations 
in rie and the random component of 1^o on scales 
much smaller than the 10.75' beam. The depo¬ 
larised canal features observed at the boundary of 
the band are consistent with being caused by RM- 
gradients at the edge of the ionised shell. 

The study presented here illustrates how even large 
objects, well-sampled by RMs, require great care to 
disentangle from confusing sources. The next genera¬ 
tion of surveys planned for the Square Kilometer Array 
and precursors instruments will enable similar studies of 
many more Galactic objects. From 2016 onwards the 






















































































Australia Square Kil ometre Array Pathfinder (ASKAP, 


Johnston et al.l [2007ll POSSUM (Polaris ation Sky Sur¬ 


vey o f the Universe’s Magnetism) project (jOaensler et alJ 
201(1) will survey the southern sky at 1 GHz and deliver 
a RM-grid with ~100 polarised sources per square de¬ 
gree 100 times the source density of the NVSS). Once 
the new catalogue becomes available it will be possible 
to identify and correct for smaller Faraday-active objects 
with greater accuracy and confidence. This work serves 
as a rehersal for future studies and highlights the chal¬ 
lenges involved. 
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A. Geometry of Simple Ionised Shell Model 

The geometry of our simple shell model is shown schematically in Fig. 1161 It consists of a spherical ionised region of 
radius R and thickness dr, containing thermal electrons at an average density rig. The edge of the shell subtends an 
angular radius (/)outer = sm~^(R/D), where D is the distance from the observer to the geometric centre. The ionised 
gas is threaded by a uniform magnetic field Bq, whose vector is confined to the y-z plane representing the Galactic 
disk. Field lines make an angle 0 to the plane of the sky in the direction of the bubble centre. From the perspective 
of an external observer, the angle C, describes the orientation of the sight-line to the y-z plane. The observed RM is 
given by Equation [3] and in the simplest case {Q = 0, Fig. [T6l- left) depends only on the line-of-sight field strength i3|| 
and the path-length L along a chord through the shell; 


R|| =-Bosin(0 - 0), (Al) 

L{(l>) = 2 y 112 - sin^ (</>)]. (A2) 

If the shell is compressing the ionised gas at the leading edge then the components of tangent to the surface will be 
amplified by a factor X. To model this effect we decompose Bq at each point on the surface into vector components 
Bzf) along the axes xi, yi and zi such that 

Bxi = Bo cos(0)sin(C), (A3) 

By,^ = Bo [ cos(0) cos(C) cos((/)i) -I- sin(0) sin((/)i) ], (A4) 

Bzi = Bo [ cos(0) cos(C) sin(0i) — sin(0) cos((/)i) ]. (A5) 

The observed magnetic field strength is then the vector sum along the line-of-sight at both the ingress {Bi) and egress 
(Be) points. From the geometry we find 

Bi = X By^ sin(a) — cos(a), (A6) 

Be = X By^ sin(a) -I- B^^ cos(a), (A7) 

where a = sin“^(0/(/)outer)- When calculating Bi the angle (fi in Equations |A3]-|A5] is sampled over the far side of the 
sphere according to <f>i = a + 4>. Similarly, when calculating Be, 4>i is sampled over the near side: 4>i = 180° — a + (p. 
The final line-of-sight magnetic field strength is then given by 

R|| = (R, + Re)/2, (A8) 

assuming Bi and Bg are constant along the path connecting the midpoint of the chord and the surface of the shell. 
The RM at any point may then be calculated by combining Equations IA2l - [A8l with Equation 01 



Fig. 16.— Geometry of an ionised shell in the near-held. We assume a constant electron density rie in the shell, which is 
threaded with a uniform magnetic held. The held is conhned to the y-z plane, representing the Galactic disk, and only the sky 
pitch-angle 0 may be varied. 



















